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LI. The Strength of Lomer—Cottrell Sessile Dislocations 


By A. N. Stro# 
Cavendish Laboratory, Cambridge + { 


[Received December 21, 1955] 


A Lomer—Cottrell sessile dislocation at the head of a piled-up group of 
dislocations may give way under the combined action of stress and 
temperature. Depending on its orientation relative to the group, the 
dislocation may yield either by recombining and slipping on a (100) plane 
or by dissociating into the two dislocations from which it was formed. The 
size of the group required for both these mechanisms is obtained as a 
function of stress, but the results depend rather sensitively on the structure 
of the dislocation core. 


§1. INTRODUCTION 
LomEr (1951) has pointed out that, in a face-centred cubic crystal, two 
dislocations in different slip planes can attract one another and combine ; 
thus the dislocation with Burgers vector $a [101] in the plane (111) can 
react with the dislocation $a [011] in (111) giving 
ga [101]+ 4a [011] > da [110]. arate. tae 


The dislocations must all be parallel to the line of intersection of the slip 
planes, i.e. to the line [110], and so the resulting dislocation $a [110] has as 
slip plane, the plane (001). Cottrell (1952) then remarked that the 
dislocation $a [110] could dissociate : 


47 (110) > 4a [112]-e 2a [110| +40 [112]; . . . (2) 
these three dislocations are all imperfect and when they separate they will 
be linked by two strips of stacking fault. This group of three dislocations 
will be unable to glide and will form an obstacle in the slip plane against 
which further dislocations can pile-up. If the original dislocations were 
dissociated into partial dislocations, 

La [101] > 4a [112]+20 [271], | 

4a [011] > $a [112]+4e [121], : 
the sessile dislocation may be formed by the direct combination of the 
second dislocation of each pair 


4a [211]-+40 [121] + 4a [110] ; Nk emeew es) 


(3) 


+ Now at Department of Physics, University of Sheffield. 
t Communicated by the Author. 
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the two remaining partial dislocations are linked to this by strips of 
_stacking fault giving the same set of three dislocations as before. We 
note that all three dislocations making up the sessile are of pure edge type. 

Other possible combinations of dislocations have been considered by 
Friedel (1955) who finds, however, that either the short range interaction is 
not such as to lead to a stable complex, or there is no long range attraction 
so that the dislocations are unlikely to meet. Thus the reaction discussed 
above appears to. be the only one of importance. 

The observation of work softening in aluminium single crystals, that is a 
decrease in the flow stress when a specimen first deformed at one tempera- 
ture is further strained at a higher temperature, has led Cottrell and 
Stokes (1955) to suggest that the sessile dislocation may give way under the 
combined action of stress and temperature. There will then be a limit to 
the size of the piled-up group of dislocations which can be formed against a 
sessile dislocation at any temperature, and raising the temperature will 
cause piled-up groups formed at the lower temperature to break up so that 
their contribution to the work hardening is lost. Friedel (1955) has also 


Higset 


\ 


—s 
(a) (b) \ 


Sessile with its apex (a) away from the pile-up (narrow stacking fault) and 
(b) towards the pile-up (broad stacking fault). 


suggested that the sessile dislocation may give way, and has interpreted 
the change from linear to parabolic hardening observed in copper single 
crystals by Blewitt et al. (1955) in terms of this. In the present paper an 
attempt is made to estimate the size of the piled-up group which the sessile 
dislocation can withstand. 


§ 2. GENERAL CONSIDERATIONS 


The piled-up group may be formed on either side of the sessile disloca- 
tion, so that the latter may have its apex pointing either away from the 
group as in fig. 1 (a), or towards the group as in fig. 1 (6). The arrows 
indicate the directions in which the stresses tend to move the partial 
dislocations forming the sessile dislocation ; thus the stacking fault will 
be narrow in the first case, wide in the second. 

If the dislocations are recombined by the reverse of reaction (2), and 
then allowed to dissociate again, the sessile dislocation may be converted 
from one arrangement to the other; the stress will clearly favour the 
dislocation’s pointing towards the group. The large activation energy 
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needed to combine the dislocations over a finite length can be avoided by 
the sequence shown in fig. 2. This represents successive positions of two of 
the three partial dislocations forming the sessile dislocation, the stacking 
fault between them being shaded. Of the dislocations, one, that with 
Burgers vector }a [110], is constrained to lie along the line of intersection 
of the two slip planes ; the other, 3a [112] or 4a [112], then crosses over it 
as shown. By repeating this process on the other slip plane the sessile 
dislocation can turn its apex towards the group. It is easily seen that the 
activation energy needed is less than twice the energy of a constriction 
W ons For dislocations of this type, the latter has been estimated by 
Stroh (1954) to be 


W 


pee ee LO ard (log dra) ete) 
where yu is the rigidity modulus and r, the radius of the dislocation core ; d 
is the distance between the dislocations and depends on the stress. This 
activation energy will be found to be small compared with that required 
for the sessile dislocation to give way; hence before this happens the 


dislocation should have turned its apex towards the piled-up group. 


Fig. 2 


(2) () (c) (d) 
In the initial uniform ribbon of stacking fault (a), a constriction is formed (6) ; 
this develops into a small loop where the dislocations have crossed 
over (c), and the length of this portion increases (d). 


When the partial dislocations cross over, the stacking fault will change 
from intrinsic (with the sequence of close packed planes abe } beabec . . . .) 
to extrinsic (abe } b abca....) or vice versa (Frank 1951, Frank and 
Nicholas 1953). These both have all nearest neighbour relationships 
correct, and the same number of incorrect second nearest neighbour 
relationships ; thus they may be expected to differ only slightly in energy. 
Since the stress due to the piled-up group will be found comparable to that 
required to extend the stacking fault by widely separating the partial 
dislocations, it should be quite great enough to overcome any small 
difference in stacking fault energies. Also any possible change in the core 
energy of the dislocations should be unimportant. On the other hand, the 
intrinsic stacking fault involves a break in the correct stacking sequence 
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across one plane, the extrinsic across two, and the possibility that the 
conversion is more difficult than energy considerations alone suggest 
cannot be entirely ruled out. In this case the sessile dislocation may give 
way from either position. 

Summing up, we may say that if the change from intrinsic to extrinsic 
stacking fault can take place with reasonable ease, the sessile dislocations 
should all turn their apices towards the piled-up group and give way from 
this position ; but if the change in stacking fault cannot occur there will be 
two types of dislocation lock, which we may expect to be of different 
strengths. 

The sessile dislocation may give way in either of two ways. First, the 
dislocations may recombine, the reverse of reaction (2), giving a perfect 
dislocation with Burgers vector 4a [110] which can then glide away on the 
(001) plane. Though this is not the slip-plane usually observed, Boas 
and Schmid (1931) have reported slip on it in aluminium at high tempera- 
tures. This mechanism will be easier when the sessile dislocation points 
away from the piled-up group, as then the partial dislocations are close 
together and their energy of recombination is less. 

Secondly, the sessile dislocation may dissociate into the two dislocations 
from which it was originally formed, which can then glide away on their 
respective slip planes. When the sessile dislocation poimts towards the 
piled-up group the dissociation can take place by the reverse of reaction 
(4); but if it should point away from the group the dislocations must 
first be recombined (by the reverse of 2) and then dissociated (the reverse 
of 1), and so the break up will be more difficult. 

Failure by recombination will be considered in § 4, and by dissociation 
in §5. In the following section, the energy of recombination per unit 
length of the dislocation is found. 


§3. ENERGY OF RECOMBINATION 


To obtain this we must first find the stresses due to the piled-up group. 
We shall need only the stresses at the position of the sessile dislocation ; 
these are readily found by a simple argument. Take y axis normal to the 
plane of the group, the (111) plane, and z axis parallel to the dislocation 
lines, the [110] direction. The Burgers vectors b of the dislocations may 
be resolved into edge components b,, and screw components b,. All the 
components of stress at a point in the slip plane due to any one dislocation 
are zero except o,., and a,, ; and these are in the ratio 


C py/Fy,=0,/(1—v) b,, A Tt. -o Ask eee Gy 


where v is Poisson’s ratio. This must also be true for the resultant stress 
due to the whole group. For the purpose of calculating the stress, suppose 
the sessile dislocation to be replaced by a dislocation similar to those of the 
group ; the group produces on it a force b, o,,+b, ¢,,. If there are 


+ This was suggested by Professor F. C. Frank. 
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n dislocations piled up under a resolved shear stress oy then according to 
Cottrell (1949), this force is nbag, or 


6. Czy tO, Fy,=NDdo5. eee «wee aay" (i) 
Solving (6) and (7) we obtain 
ae b,,bnoy 
ey ie Ae eel ee eel) 


With b= $a [101]=a [112]+ 4a [110], eqns. (8) reduce to 
Oxy=2V/3no,/(4—v), o,,=2(1—v)no,/(4—v). . . . (9) 


We also need the stresses on the (111) plane ; using dashes to denote 
corresponding quantities for this plane, we obtain by the usual trans- 


formation 
# 


eee! a= 
—9 Fay Oyp=3 7 (10) 


The stresses (9) lead to a force per unit length on the Wey 
da [112] of 1/2 no, a/(4—v), while (10) gives a force 7\/2no,a/9(4—v) on 
4a [112]. If the widths of the two ribbons of stacking fault are d(1—€&) 
and d(1+-&), so that d is the mean width, then the energy per unit length 
due to the stress and to the stacking fault energy y is 


2d 4164/2 0 o, ad (I—46)/9 (4-2), . . . . (11) 


where the positive or negative sign must be taken according as the stress 
tends to decrease or increase the width of the stacking fault, that is 
according as the apex of the sessile dislocation points away from or towards 
the piled-up group. To (11) we must add the energy of interaction of the 
dislocations. 

Now according to Nabarro (1952) the energy per unit length of two edge 
dislocations with Burgers vectors b’ and b” separated by a vector r is 


ere toe D0) (Decr) 


Summing contributions of this form, we find the energy of interaction of 
the three partial dislocations of the sessile dislocation is, 


pa d Nee As } 
= ty es + 2 log Gap A re yrieg caf (13) 


where a constant has been added so that the energy is zero when the 
dislocations are recombined. If <1, (13) becomes 


pa 2 4 se : 14 
mes =a) {logs +42--O(E . (14) 
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The total energy is the sum of (11) and (14) ; differentiating this sum 
with respect to d, we find the energy to be a minimum when ‘Lo 


De gid, 5. 169 pee i oh ye 
a 360 eae 2) 8 ane 
The corresponding value of the energy is 
y LN MOSEL Fy aaa 1) ee eae 
AP st (loge + é 4° oS 


We have now to vary é so as to minimize (16). When the upper sign in 
(15) is taken we find that for all values of nop, €<1/64; thus not only is 
our original assumption that €<1 justified, but we may also to a good 
approximation neglect all terms in €. Terms in may also be neglected, 
though to a less good approximation, when the lower sign in (14) is taken 
provided no a is small compared with y. Since the stacking fault on (111) 
becomes infinitely wide (or rather the dislocation 4a [112] is displaced a 
large distance from the piled-up group to where the stress is less than the 


value (9)), when . 4/2n6,a/(4--v) 95. a, ee ee 


the condition implies that the dislocations must not be too widely separated 
by the stress of the group if it is not fulfilled the energy of recombination 
becomes so large that the process is no longer of interest. Neglecting &, 
we may write the energy of recombination 


T= {pa7/36n (l—v)} Ly, oo. a 
where | L,=log (d/r))—1. 


When the stress is so large that d is comparable to the radius of the core 
1, this expression cannot hold, as 7” must always be positive if the extended 
state is stable. The difficulty arises because the interaction between the 
dislocations is no longer given correctly by elastic theory at short distances. 
If, to ensure that the force between two dislocations remains finite, we 
assume it varies as (7?-+-¢?)~1/2 where @ is a constant length (the width of 
the dislocation), instead of as r~1, then in the expression for the energy 
(12), we must replace log (r/rg) by log {r+-(7?+@)12V/¢. With this the 
width of the stacking fault is given by 


(d?+-€?)¥/2— {ua?/727 (1—v)} (y+ 0-348n of a), 


and the energy of recombination per unit length is still given by (18), but 
where now 
L,=log {d+(2+2)!2}/6—dl(@+2)12 ; 
if 7 is defined by 
cosh n= (d2+ £2)1/2/¢— ONCE 


(y/pa)+0°343 (no9/1)’ 
then Lye tanh yee ir See a ee oe OY 
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The stacking fault energy y has according to Seeger and Schoeck (1953) 
the values 200 erg/cm? for aluminium, and 40 erg/em2 for copper ; the 
dimensionless quantity y/.a then becomes 1-86 x 10-2 for aluminium, and 
0-24 x 10 for copper. On account of the high value for aluminium, the 
width of the stacking fault will be narrow ; if the dislocation width C is 
greater than 0-35a, (19) indicates that the dislocation will be completely 
recombined even in the absence of any applied stress. As the dislocation 
will not then be sessile, we must suppose £<0-35a. On the other hand it 
seems unlikely that ¢ will be much smaller than the Burgers vector which is 
a/,\/6=0-4la. These considerations suggest that 0-2a or 0-3a are suitable 
choices for ¢. For copper larger values of € would still lead to an extended 
dislocation ; nevertheless in numerical calculations the same values will 
be used so that the results may more readily be compared. 


Fig. 3 


XN vA 
XN / 
‘\ 4 
N / 

N 4 

4 
EES / 
2. ‘\ 4 
Nec, 
20 
‘\ 7 
v 
e) 


A dislocation recombined over a length /, and bowing out on a (100) plane. 


§ 4. FarnurE BY RECOMBINATION 


Suppose that the dislocations are recombined over a length /; this 
length of dislocation can then bow out on its slip plane, the (100) plane, in 
a manner reminiscent of a Frank—Read source. If o’ is the resolved shear 
stress on this plane, the radius of curvature Ff is given by 

Rone xy Seely Te eee Se pyc ts, ws (2) 
where a/,/2 is the Burgers vector and 7’ the line tension of the recombined 
dislocation. If the dislocation subtends an angle 20 at its centre of curva- 
ture, the energy required to bow it out is made up of the following 
contributions : 

(i) Twice the energy of a constriction, the factor 2 occurring since 
W constr ftom (5) is the energy of a constriction in only one of the stacking 
fault ribbons ; 

(ii) The energy of recombination of a length /, T'/; ; 

(iii) The energy due to increase in length of the dislocation line, 
T (2R6—1) ; 
and (iv) minus the work done by the stress (a/\/2)o° (R?6—4R? sin 26). 
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Hence the total energy is 
W=2Weonste + LI+T (2RO—1)—(a/1/2) oR? (8@—} sin 20) 
= 2W +27’R sin O6+T7R (@—2sin6+}sin 26), . ... ; (22) 


constr 
on using (21) and /=2R sin 6. 

W, is a maximum for variations of 6 when 
cos 0=1—T"/T ; eS ae 


this gives a configuration of unstable equilibrium, and from this position 
the loop can grow without requiring further help from thermal fluctuations. 
Now the line tension 7’ is 


= {ua?/8m (1—v)} log (R/ro), - - - - + = (24) 


and comparing this with (18) we see 7”/7'<1; then (23) gives approxi- 
mately 0=(27"/T)/?. Substituting this in (22) we obtain the activation 


energy W1=2 Ween +8122 TUY30'a. . . . . . (25) 
From the stresses (9) we obtain o’=40,,,=2noy/ 1/3 (4—v) ; then (18), (24) 
and (25) give 


7, a 243 Ly 3/2 i 1/2 
Wi=2 at W ponseale \ vt 


1084/6 7? (1—v)? no 


=2W sonete + 3°16 X 10-2 wa? (u/noy) Ly LeM?, . . (26) 
where Le=log Rjr,~loe (uinanye - seals ee 


The probability of the sessile dislocation’s failing in time ¢ is 
vt (ja) exp (—W,/k7), 


where v is an atomic frequency of vibration, and / the length of the 
dislocation, so that //a is the number of points at which the breakaway can 
start. Equating this probability to unity we find the dislocation will give 
way when 
Wi=kT log (vil/a)s a a eed 


with v~1012 sec —1, t~108 sec, l/a~104, this is W,=44 kT. 

From (26) and (28) the stress at which the sessile dislocation gives way 
may be found as a function of temperature ; results are shown in fig. 4, for 
the case in which the dislocation points away from the piled-up group. 
As typical results we may note that under a stress of 1 kg/mm? the 
maximum number of dislocations in the group varies for copper from about 
150 at 300°K to 220 at 100°K ; for aluminium at the same two temperatures 
the numbers are 25 and 45. 

When the dislocation points towards the group, the activation energy 
will be large for large stresses because the stacking fault is wide, and also 
for small stresses because the size of the critical loop is then large ; thus 
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there will be a minimum activation energy for which the dislocation can 
give way. For aluminium this minimum is 2:6 ev if £=0-2a and 0:80 ey if 
¢=0-3a, the corresponding temperatures according to (28) being 680°K and 
210°K ; the possibility of the sessile dislocation’s failing in this way thus 
depends very sensitively on the structure of the dislocation core. For 
copper, the greater width of the stacking fault leads to so great a value for 


the energy (of the order of 100 ev) that this mechanism can be altogether 
ruled out. 


Fig. 4 


0:04 | 


0-03 


NG of L-> 


0-02 


0-01 


Stress at which the sessile gives way by recombination plotted as a function 
of temperature for two different values of the dislocation width ¢: 
(a) copper, C/a=0-2, (b) copper, ¢€/a=0-3, (c) aluminium, ¢/a—0-2, 
(d) aluminium ¢/a=0:3. 


§5. Faure By DissocraTion 


As was mentioned in § 2, the sessile dislocation can dissociate more easily 
when it points towards the piled-up group ; we shall suppose this to be the 
case. If the number of dislocations in the group is greater than indicated 
by eqn. (17), the dislocations 4a [112] and 4a [112] will be repelled to large 
distances, and so will not influence the dissociation of the dislocation 
da [1 10]. When this dissociates, by the reverse of reaction (4), work must 
be done in separating the two dislocations so formed, and the separation 
can take place most easily by first forming a small loop (fig. 5). Since 
forming such a loop involves doing work against the line tension, we expect 
the energy to be least when only one of the two dislocations deforms to any 
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great extent. The stress on }a [211] is the greater, and so this is the 
dislocation which can deform more easily ; the dislocation ja [121] is 
assumed to remain straight and undeformed along the line [110], the z axis. 
The energy of the loop can be found by a method similar to that used by 
Stroh (1954) to obtain the energy of a constriction ; this gave values of the 
constriction energy in good agreement with those obtained by Seeger and 
Schoeck (1955) by a variation method. 

A small element of the dislocation 4a [211] of length 6s at the point 
P (x, z) is subjected to the following forces : 


(i) A force due to the line tension 7' (d?x/dz*)és where 


par 1 R 
== a — > 8 ee 
{1+ saan lees, ; fo 
(ii) A force due to the attraction of the dislocation 4a [121] which we 
take to be — Ads/(x?+ C?)1/2, where 


A=(pa4/167){1+2 (I—v)“4};  . . . . ~ (80) 
and (iii) a force F'ds due to the stress and to the stacking fault where 
F=y-+no,a (2—v)/+/2(4—), Meer pres (CRY 


and we have used the stresses (9). 

The energy of the loop will be greatest when it is in a position of unstable 
equilibrium ; for equilibrium the above forces give 

T (d*2x/dz*)— A(x?+ (?)-V2+4+ F=0. 
Integrating, we have 
47 (dx/dz)?=A log {a/€+(1+a?/C?)/2}—Fa, . . . (82) 

since da/dz=0 when x=0. At the pomt z=0 at which the separation 
is greatest dx/dz is zero, so that the maximum separation x,, can be found 
by equating the right hand side of (32) to zero and solving for x ; the result 
is plotted in fig. (6) as a function of FZ/A. 

Now we have the following contributions to the energy of the loop : 

(i) $7’ f (dx/dz)? dz due to the line tension ; 


(ii) A flog {a/{+(1+-a?/¢?)1/?} dz due to the mutual attraction of the 
dislocations ; 


and (iii) —F’ f x dz due to the stress ; 


where in each case the integration extends over the whole length of 
dislocation which has dissociated. The total energy is then 


W.=J [$7 (da/dz)?+-A log {a/€+-(1-+-a?/L2)1/2}— Fa] dz 
=2 (27A)WV2 in [log {a/f+-(1-+-a?/f?)1/2}— Fy] A]1/2 da 
on using (32). Substituting x=¢ sinh wu we obtain 
Wy=2 (27 A)UIT (PCA he ame ges 
where oS be (w—v sinh w)!/2 cosh wu du 


1 


Um 
=> | (w—v sinh u)1/2 du, er eee (ee § 
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and w,, is thé positive value of w for which the integrand is zero; the 


function vf (v) is plotted in fig. 7. Equations (29), (30). (31 
with v=1/3, q (29), (30), (31) and (33) give, 


W.=5-26 x 10? aE (log a,,/C)1/2 f CRAY 5 ete (35) 
and PC/A= {48-0 (y/ua)-+13-8 (noptei(Gid).. t. 2. . (36) 
Fig. 5 
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% @ [110] dislocation dissociated over a short length. 


Fig. 6 
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Size of critical loop as function of stress. 


The above analysis holds provided the size of the critical loop is small 
compared with the width of the stacking fault. When numerical values 
are inserted, this is found to be so for copper for all stresses giving an 
activation energy of reasonable magnitude; in this case we may deter- 
mine the stress at which the sessile dislocation breaks down as a function 
of temperature from eqns. (28) and (35). For a metal with a high stacking 
fault energy such as aluminium this will only hold at high stresses (low 
temperatures) ; for lower stresses the size of the critical loop will be greater 
than the width of the stacking fault. If it is much greater than the 
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stacking fault width, we can neglect the dissociation into partial dis- 
locations, and suppose the dislocation breaks up by the reverse of reaction 
(1). Then we may use the previous treatment and the energy is still 
given by eqn. (33) but where now 


A=(pa/167){1-++1/(1—»)} , 

T = (y0?/327){1+3/(1—v) } log (@,,/¢’), 
and F=no, a/v/2; . ae 
we do not assume ¢ has the same value as before and so denote it now by ¢’. 
With v=4 we obtain the energy 


(37) 


W,=0-15 pa? C’ (log x,,/6')1/? f (28-40 op Cjua). - . ~ (38) 
Fig. 7 
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Value of the function ofr) =|" (u—v sinh wu)!” du. 


When the size of the critical loop and the width of the stacking fault 
are approximately equal neither eqns. (36) nor (38) apply. However, 
since both the size of the loop and the width of the stacking fault vary 
rapidly with stress at this point, the range of stress for which neither 
equation is valid is small and may be ignored. The curves for aluminium 
in fig. 8 have been obtained by using eqn. (36) when it is less than the width 
of the stacking fault and (38) when it is greater, and choosing ¢' so that 
the energy varies continuously. Under a stress of 1 kg/mm? the maximum 
number of dislocations in a piled-up group varies for copper from about 
500 at 300°K to 750 at 100°K ; for aluminium at the same temperatures 
the numbers are 200 and 350. These numbers are several times those 
obtained in §4, and illustrate the differences in strength of the sessile 
dislocation in the two situations, 
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When the sessile dislocation points away from the piled-up group, the 
activation energy for dissociation will be greater than if it pointed towards 
the group. Then comparing figs. 4 and 8 we see that the dislocation will 
give way by recombination before dissociation can occur. Accordingly 
there is no need to consider dissociation in detail in this case. 


Fig. 8 


0 100 200 300 400 500 600 700°x 
T= 
Stress at which the sessile gives way by dissociation plotted as a function 
of temperature for the different values of the dislocation width ¢: 
(a) copper, C/a=—0-2; (b) copper, ¢/a=0-3; (c) aluminium, ¢/a=—0-2; 
(d) aluminium, ¢/a—0-3. 


§6. Discussion 


As we have seen the results depend somewhat sensitively on the structure 
of the dislocation core, and accordingly not too much weight can be 
attached to the exact numerical values obtained. We can, however, draw 
some general conclusions. When the sessile dislocation points away from 
the piled-up group it will give way most easily by recombination and slip 
on (100); in this position it provides a rather weak obstacle, and in 
aluminium especially only a very small group would be needed to raise the 
stress sufficiently to overcome the sessile dislocation. If the dislocation 
can change its orientation it will strengthen itself considerably by pointing 
towards the group. If it cannot so change we will have two types of 
obstacles of rather different strengths; sessile dislocations pointing towards ° 
the piled-up group, being the stronger, will be the more important. From 
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this position the dislocation will give way by dissociation, except at high 
temperatures in metals with high stacking fault energy, such as aluminium, 
when recombination will be easier. At the temperature at which the 
latter begins there will be a marked drop in the size of the pile-up which 
can form. 

Finally, we note that in no case can particularly large piled-up groups _ 
be formed. For comparison we may use the result of Stroh (1954 b) 
that to initiate a crack we require a piled-up group of magnitude given by 
no,/u=1-5, and this is greater by a factor of 10 or more than the groups 
which the sessile can withstand. We suggest that the reason why single 
crystals of face-centred cubic metals never fracture but always draw down 
to a point is that there are no obstacles which will withstand a sufficiently 
great piled-up group ; in polycrystals, on the other hand, grain boundaries 
can act as the barriers. 


,s 
ACKNOWLEDGMENTS 


It is a pleasure to thank Dr. A. H. Cottrell for suggesting the problem, 
and Dr. J. Friedel for a number of stimulating discussions. The work was 
done during the tenure of a D.S.I.R. senior award. 


REFERENCES 


Biewitt, T. H., Cottman, R. R., and Repman, J. K., 1955, Defects in Crystalline 
Solids (London: The Physical Society), p. 369. 

Boas, W., and Scumip, E., 1931, Z. Physik, 71, 703. 

CorTRELL, A. H., 1949, Progress in Metal Physics 1; 1952, Phil. Mag., 48, 645. 

CorTrELL, A. H., and Stroxss, R. J., 1955, Proc. Roy. Soc. A, 233, 17. 

Frank, F. C., 1951, Phil. Mag., 48, 1224. 

Frank, F. C., and Nicuouas, J. F., 1953, Phil. Mag., 44, 1213. 

FRIEDEL, J., 1955, Phil. Mag., 46, 1169. 

Lomer, W. M., 1951, Phil. Mag., 42, 1327. 

Nasarro, F. R. N., 1952, Advances in Physics, 1, 271. 

ScHorck, G., and Sercer, A., 1955, Defects in Crystalline Solids (London : 
The Physical Society), p. 340. 

SEEGER, A., and ScHorck, G., 1953, Acta Met., 1, 519. 

Srrou, A. N., 1954a, Proc. Phys. Soc. B, 67, 427; 1954b, Proc. Roy. Soc. A, 
2238, 404. 


LI. The Long B-Lifetime of 4C and the 4*N Spectrum 
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ABSTRACT 


The long -lifetime of C and the spectrum of !N are explained using 
a small tensor force in addition to the usual central and spin-orbit forces. 
The y-decay branching ratio of the 3-95 Mev level of 14N is predicted and 
also the nucleon reduced widths for the low levels of 4 and 4N. 


§ 1. LyrRopuctTIoNn 


THE intermediate coupling shell model has met with such success (Inglis 
1953, Zeldes 1953, Lane 1953, Elliott and Flowers 1955) in explaining 
the positions and properties of low levels in light nuclei that one welcomes 
any violent disagreement between theory and experiment as a means of 
learning more details of the model. The B-decay of 4C to 14N is such an 
apparent disagreement. According to the selection rules and the shape 
of the spectrum the transition is an allowed one, but the decay constant is 
such that log ft=9-03. Thus we have a transition which, although 
formally allowed, is hindered by a factor of the order of 10*. 

Bouchez (1950) suggested that the large ft value results from a change 
in the orbital angular momentum which breaks the selection rules for an 
allowed transition (L-forbiddeness) but, on its own, this would demand a 
degree of L-S coupling quite incompatible with that of neighbouring 
nuclei. The explanation by Inglis (1953) depending on strong admixtures 
of higher configurations is also hard to accept in view of the overall success 
of the theory based on the lowest configuration. 

Inglis (1953) showed that no mixture of central and spin-orbit forces 
could produce the cancellation in the matrix element necessary to give 
such a large ft value, but his remark that the introduction of a tensor force 
also failed to give cancellation was incorrect. Jancovici and Talmi 
(1954) showed that the inclusion of a tensor force with the central and 
spin-orbit forces of intermediate coupling gave the desired cancellation 
but that one needed such a large tensor force that the resulting spectrum 
of 4N was in violent disagreement with experiment. In particular it 
gave a ground state with total angular momentum /=2 instead of the 
observed J=1. 


wee A I ee ee Ee 


+ Communicated by B. H. Flowers. 
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The purpose of this paper is to show that the trouble which Jancovici 
and Talmi experienced with the spectrum arose from the ratio of two 
tensor force radial integrals which they calculated. This ratio turns out 
to be very sensitive to the radial dependence of the interaction and of the 
wave functions, both of which are uncertain. We treat this ratio as a 
parameter and show later that the value necessary to give the observed 
spectrum is a very reasonable one. Having determined wave functions 
which give the correct spectrum we go on to calculate relative reduced 
widths for the reactions “N(d, p)&N and “C(d, n)4N and the branching 
ratio for the y-decay of the second excited state in MN at 3-95 Mev. 
If these also agree with experiment our procedure is largely vindicated. 


§2. Tor InrERAcTION Matrix ELEMENTS 


A general two-body interaction consisting of a central force together 
with spin-orbit and tensor terms may be written 
Mie Bat pet ibS Netter Ube Era I;) 
<j t 
+ ET (ij) (toler) ~ oe) Vira). 
<j Un) 3 
The matrices of this interaction in the configuration p!° appropriate to 
mass 14 may be worked out by the standard processes described, for 
example, by Elliott (1953). They are given below for each value of J 
and of 7’, the isotopic spin :— 
11) 
138 wd DE 3 ad eet 


(L+2K) 3/5 C Vie 
(L—K)+4é—-214 | —vVé 
(L—3k) s 


T=1 
ata 3p, 


T=0 | (L-+2K)q V2é 
(L—3K)t+£+10B 


eal BP 33P 


J=2 |) (L—K)+4 £4214 of =e) (L—3K)t+4é—5B 


SAD® 317), 33P, 


J=3 | (L—K)—é—64 J=2 | (L—K)g —/xé 


(L—3K)t—}E+B 
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Here L and K are the usual central force radial integrals for the p shell, 
A, B and C are similar tensor force integrals which we shall discuss later, 
and q, s and ¢ describe the exchange character of the central force, viz. :— 


qg=W+M—H—B, s=W—M+H-—B, t=W—M—H-+B, with 
W+M+H+B=1. 


The factor 7'(ij) in (1) describes the isotopic spin dependence of the tensor 
force and is of course contained in the parameters 4, B and C. We 
write [«, 6, y] for the vector corresponding to the lowest root of the 
J=1, T=0 matrix, and similarly [u, v] and [6, 4] for the (0, 1) and (2, 1) 
matrices respectively.f 


§3. THE 6-Decay Matrix ELEMENTS 


Since the “C(8~)4N transition involves a change of spin the Fermi 
matrix element vanishes. Calling the Gamow-Teller matrix element G, 
one finds 


GFONAEAOVE Tees oe 1. 3 ot (2) 


In order to obtain log ft9, this expression must be so small that we can 
assume it to vanish. Now, to produce the right sign of the doublet splitting 
in 5N we must have €<0; hence from the (0, 1) matrix it follows that 
(w/v) >0 for the lowest root. A necessary condition for the vanishing of 
G is therefore that («/y)<0; and precisely, 


(ES VES are a ee aC) 


If C=O, as in the absence of a tensor force, it can be seen from the (1, 0) 
matrix that («/)>0 for the lowest root and so G cannot vanish. Further, 
it can be seen that the vanishing of G depends sensitively on the value of C. 


§ 4. GENERAL RESTRICTIONS ON THE PARAMETERS 


The general matrices written down in § 2 contain a large number of 
parameters but there are two at least that can be estimated with some 
certainty. They are the ratios (L/K)=6 and (€/K)=4. The first is a 
slowly varying function of range and potential shape ; the latter follows 
both from the intermediate coupling work of Lane (1953) on the neigh- 
bouring nucleus #C on the one side, and from the doublet splitting of 
15N on the other. 

If we now write 


80=—27K" and ~ 2fA=yrk 
we can treat K as an energy unit, leaving only three parameters, x, 7 


and s in the 7=0 matrices. Instead of s we find it convenient to use 
the parameter 4, being the excitation in units of K of the (2, 0) level 


ae ee eee ee 


+ We use the notation (J, 7’) to indicate spin and isotopic spin of a level. 
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above the lowest (1, 0) level. This has the added advantage that we can set 
a lower limit 4K >6 from the observed spectrum. Let us now consider the 
(1, 0) matrix, making the condition that («/y)<0. This is a necessary, 
but insufficient, condition that the matrix element G should vanish. 
For each value of y we can then find regions of the (x, 4) plane which satisfy 
this condition by drawing the curves of <=0 and of y=0. The regions 
are shown cross-hatched in figs. 1, 2 and 3 for values of y of 1, 2 and 4 
respectively. Furthermore it can be shown that the double-hatched 
region labelled A is the only one allowed because that is the only region in 
which the negative value of (z/y) corresponds to the lowest root of the 
matrix. It must be remembered that for each point x, 4 the exchange 
parameter s is determined and so there will be a further restriction if we 
have a priort knowledge of s. 
Fig. 1 


Allowed region A of (x, 4) plane for y=1. 


The figures show that « must be positive, which gives a tensor force of 
the same sign as that needed for the deuteron. For small y there is a 
large range of x for which4 is reasonable in the allowed region, while for 
larger y one has to come to very small x to get a large enough J. 

Such a small value of x is hardly acceptable since it entails a large nega- 
tive value of s, incompatible with the usual exchange mixtures suggested 
by Rosenfeld (s=—1-8) and Serber (s=0), It follows that for a large 


Fig. 2 


FF 


) . 


Allowed region A of (x, 4) plane for y=2. 
Fig. 3 


(\ A 


x<=O0 


S=0 
Allowed region A of (x, 4) plane for y=4. 
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enough J and reasonable s we must look to the small values of y. Jancovici 
and Talmi (1953) used a large value of y with a reasonable s resulting in a 
low (even negative) value of 4, and this gave the (2, 0) level as ground state 
in disagreement with observation. 


§5. ASSUMPTION OF CHARGE SYMMETRIC INTERACTION, LEADING TO 
SPECTRUM DETERMINATION 


In order to insert the precise condition (3) for the vanishing of G one 
must solve the (0, 1) matrix as well as the (1, 0) matrix, involving three 
new parameters g, f and B. It is fortunate that they can be related to 
those parameters already considered if we introduce the single assumption 
of charge symmetric forces, i.e., M=2B, H=2W and T(ij)=(7;.7;). For 
the remainder of the paper we make this assumption, whereupon 


g=—3s,t=—4 and 3B=C—7A. 


If we include the precise condition (3) we are able to plot the spectrum in 
units of K against 2, the measure of strength of the tensor force ; and this 
ean be done for each value of y we consider. Figures 4, 5 and 6 show all 
the calculated levels of the configuration p1!° in the first 64 of energy 
excitation, which are to be compared with the observed spectrum in 
fig. 7, from Ajzenberg and Lauritsen (1955). 

Above 4 Mev the observed spectrum is uncertain but we can be 
reasonably sure of the properties of the ground and lowest two excited 
levels. The calculated levels fit it very well with this pattern if we take 
x close to 0-4 in each figure. Fitting the first two excited levels determines 
x and K for each y, and g and 4K then follow in table 1. 


Table 1. Parameter Values Deduced from Fitting Spectrum 


1 5 0 
Fig. 5 2 0-395 —1:8 0:70 5-8 
Fig. 6 4 0-384 —2-5 0-70 45 


The calculated excitation 4h of the (2, 0) level is seen to decrease as y 
increases, and since no (2, 0) level is observed in the first 5 Mev we must 
keep y small. It is probable that the (2, 0) level lies in the group observed 
between 4:9 and 6:4 Mev which would be consistent with y=1 or 2. 
The (2, 1) level, probably the analogue of the 6-7 Mev level in ™C, lies in a 
reasonable position around 10 Mev. The scarcity of calculated levels 
near 6 Mev implies that all but one of the even parity levels observed there 


ee SS Ee 
+ In the alternative notation for the central force, @g,/4,=(1+39)/(3 —3q). 


Calculated 14N spectrum relative to lowest (1, 0) state as zero, with y=l. 


Fig. 5 
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Calculated 14N spectrum relative to lowest (1, 0) state as zero, with y=2 
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Calculated *‘N spectrum relative to lowest (1, 0) state as zero, with y=4. 


3-95 J=l, T=O 
2:31 sic T=! 
O J=1l, T=O 


Observed !4N spectrum. 
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must come from some higher configurations, but this is to be expected at 
such an excitation. The odd parity levels cannot of course come from 
the lowest configuration p!° considered here. 


§ 6. Discusston or DepucED PARAMETER VALUES 


We must now ask whether the parameter values of table 1 are reasonable. 
The value g=0-7 which comes out for each value of y is in good agreement 
with the value g=0-6 of the Rosenfeld central force based on the deuteron 
singlet-to-triplet strength ratio. It is found to be so constant here because 
we have fitted the (0, 1) level excitation above the ground state (1, 0), and 
this is sensitive to g. In his review, Inglis (1953) favoured a value of 


Fig. 8 


(iii) 


- (*/) 
O 0:5 1-O is 2:0 
The ratio y=—7A/C for various interactions and oscillator wave functions 
ocr e777/2b%, 
(i) Yukawa e~/¢/(r/a). 
(ii) Gauss e-""/@, 


ae ae 1 eee a ,@ 
(iii) Field theoretic —— Ee +3 = | : 
ria i i 


K=-—1-3 mev which is in agreement with that found here for small y 
values. It is not simple to discuss the values of x and y because they 
both depend on the integral C' which we find to be very sensitive to our 
choice of the radial dependence of the interaction and wave functions. 
With a Yukawa interaction of the usual range, and oscillator wave 
functions, C comes out as a difference of two nearly equal integrals and 
is therefore small. Thus in order to get a sufficiently large C with such 
assumptions one must take a very strong tensor force, making the integral 
A, which does not involve a cancellation, very large.. This disrupts the 
spectrum as in the work of Jancovici and Talmi (1953). The key to the 
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problem therefore appears to lie in a larger value of (C/A), i.e. ina smaller 
y. The ratio 
—7A 
fs tEe 
is plotted against a/b in fig. 8 for various interaction shapes, where a 
is the range of interaction and b the wave function parameter. Although 
a and bare uncertain we can make rough estimates of a=1-37 and 
b=1-64 (in units of 10-13 cm) from the 7-meson mass and nuclear size 
respectively. 

It is seen that although the Yukawa force gives a large y for reasonable 
a/b values, both the Gauss and the singular field theoretic forms give 
smaller y. However, y varies so much with a and 6 that one can say 
little more than that a small value of y such as is needed here is not 
unreasonable. 

The strength of the tensor force is best illustrated by saying that the 
largest diagonal matrix element in this problem is 214=—0-6 Mev 
for y=1, while the coupling term is 3\/5 C=1-3 mev. The spin-orbit 
coupling terms are of the order of €=6 Mev, while the central force 
differences between diagonal terms is of the order of 3Kh—4-:5 mev. | It 
can therefore be said that the tensor force in this problem is only a 
perturbation on the intermediate coupling result. The reason why a 
perturbation has such a large effect on the f-lifetime for this nucleus is 
discussed in § 10. 

Although, as we have said above, it is difficult to relate the values of 
these parameters to the interaction constants because of the arbitrariness 
of the radial shapes, we can nevertheless obtain rough estimates. Let us 
assume the functions 

An eee nD leet 2) exp (—r;,/4) 
V rid=V, ot Haseee » Virg= y ea 
with oscillator wave functions u(r) ocr exp (—r?/2b?), and set a=b. Since 
C is so sensitive, we use A and K to obtain the estimate ; and taking 
the parameter values from the y=1 set of table 1 we find 


V ~—60 Mev, V ,~20 Mev. 


We have already commented that the tensor force needed here has the 
same sign as that needed for the deuteron. The strength may be com- 
pared with the figures given by Pease and Feshbach (1952) who fitted the 
low energy, two-body data and the triton binding energy. For their best 
fit they found 


V ==—47 Mev, V;—24Mev; 


and although these are based on ranges different from those used here 
the general agreement in sign and magnitude is encouraging. 

Our value of V’, is proportional to y and depends sensitively on the 
position of the (2, 0) level which has not yet been assigned experimentally. 
It is therefore of great interest to know where this level lies. 
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§ 7. THE MAGNETIC anD QUADRUPOLE MomENtTS 


Having determined the wave functions for the (1, 0) ground state in 
§ 5 it is a simple matter to calculate the magnetic and quadrupole moments 
(see for example Elliott and Flowers 1955). They are found to vary very 
slowly with x and y, and using the y—1 set of parameter values from table | 
the results are compared with experiment in table 2. The agreement is 


Table 2. Magnetic and Quadrupole Moments of !4N 


Calculated | Observed 


p(n.) 0-32 0-40 
(10-6 em?) 1-2 Ibstl 


as good as is meaningful in view of the uncertainty in these quantities 
which arise from exchange currents, etc. It is pertinent to point out 
that because the magnetic moment is so insensitive to changes in the 
parameters an approach such as that of Sherr et al. (1955) which determines 
the wave function from an exact fit to the magnetic moment is open to 
large errors. A small correction to the magnetic moment due to exchange 
moments or collective effects could produce a large change in the deduced 
wave functions. 


§ 8. THE y-DECAY OF THE 3-95 AND 2:31 Mev LEVELS oF 4N 


The (1, 0) level at 3-95 Mev (see fig. 7) can decay by an M1 transition 
either to the ground state (1, 0) or to the (0, 1) level at 2-31 Mev. 
In addition it can decay to the ground state by an E2 transition. 
Let us first consider therefore the branching ratio for pure M1 radiation 

P(1*1) - 
Awa T(1*50) 

where I’ is the width of the level for such a transition. From energy 
arguments alone one might expect Ay,>1. However, it is a simple 
matter (cf. Elliott and Flowers 1955) to calculate A,,, correctly from the 
wave functions deduced in § 5, with the striking result that Ay, -<0-04 for 
all values of 2 and y considered. Those values for y=1 in table | give 
dyn =0-003. Although for low energies an M1 process is generally more 
probable than an E2, the inhibition of the 1*-1 transition is so great 
that the two processes are comparable. In fact, calculating the E2 
width for the 1*+1 transition in a similar way we find it to be twice 
as great as the M1 width. The combined branching ratio is then 


Avene & 9-01 


where the 1*+1 transition is predominantly E2. <Any collective motion, 
not considered here, would enhance this E2 transition and so increase A. 
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Experimentally the situation is uncertain. Thomas and Lauritsen 
(1952) using the 18C(d, n)!4N reaction found no definite 3-95 Mev gamma- 
ray. Woodbury, Day and Tollestrop (1953) using *C(p, y)4N estimated 
\=0-5, but the complex nature of the spectrum makes this rather 
uncertain. Clegg and Wilkinson (1953) from the same reaction did not 
see the 3-95 Mev gamma-ray and set an upper limit of rather less than 0-5. 
Mackin, Mims and Mills (1955) using C(d,n)!4N assigned a gamma-ray 
to this transition with a branching ratio consistent with Woodbury et al. 
Clearly it would be of interest to have some more definite information on 
this decay which would then throw light on collective effects. 

Sherr, Gerhart, Horie and Hornyak (1955) have given an upper limit 
of 3:5 10-13 sec for the half-life of the (0, 1) level at 2:31 Mev in 14N. 
For the y=1 set of parameters from table 1 we find a half-life of 
2-6 10714 see which is consistent with their limit. Again, an actual 
determination of this lifetime would be valuable. 


$9. THe Nucteon RepucED WiptTHs oF Low LEVELS 
or 44N anp )N 

Consider the nucleon reduced width for the break up of a state Y into 
a state Y of one less particle and a nucleon of orbital angular momentum 1. 
For simplicity in writing, and since it will be sufficient for our present 
purpose, we assume that both ¥ and Y belong to definite configurations. 
Then the width is non-zero only if the configuration of ¥ is that of ¥ 
with one extra /-particle. Suppose that in ¥ there are a number of closed 
shells and » equivalent /-particles. Then. in the notation of Lane and 
Thomas (1955), the nucleon reduced width, in units of the single particle 
width is 


At=n | Coc, PI +VZ(2x-+1) | 2 {(28-+1)(2L-+ 1}? 


vy 


x CHW) Gl dD) Cb) 5 , | {° B ale co ee 


Here (|) are expansion coefficients of the initial state in terms of 
the complete L-S coupled set %; <%{|%s) are the fractional parentage 
coefficients for the reduction /"+/"-1+-1; and dip | WP) are expansion 
coefficients for the final state in terms of the complete L-S coupled set ub. 
The + sign refers to neutron or proton respectively which are treated as 
equivalent particles in the isotopic spin formalism. The Wigner 6-j 
symbol is defined, for example, by Elliott and Flowers (1955). Although 
the absolute values of widths extracted from stripping theory are not 
reliable, the relative widths of various levels of the final nucleus should be 
meaningful. 

If we interpret the ground and 6-3 Mev levels of !®N as the | /2- and 3/2— 
doublet of the configuration p!! the wave functions are unique. In 
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table 3 we give the ratio of the reduced width of these two levels for 
break up to the ground state of 4N using (4). The first three columns 
give values for the ratio which follow from 4N wave functions corre- 
sponding to extreme coupling models, while the fourth gives the value 
resulting from the parameters of the y=1 set in table 1. 

This value is not sensitive to changes in y and corr esponds to an absolute 
value of A?(1/2-)=1-4. Experimentally no figure is available although 
it could be found from the reaction N(d, p)®N. 


Table 3. A2(3)/ A2(4) for 5N with Various Wave Functions 


ies ee L-S ae 5-j extreme ipl ears 
1-00 0-10 0 0-05 


If we take the °C ground state wave function from the work of Lane 
(1953) we can similarly calculate the relative widths of various levels in 
MN for the reaction C(d, n)4N. The widths relative to the ground 
state are given in table 4 together with the extreme coupling values. The 
absolute value of A?(gs) corresponding to column 7 is 0-68. To our 
knowledge, experimental results are not yet available for these widths. 


Table 4. A?/ A? (Ground State) for 4N with Various Wave Functions 


L-S extreme j-j extreme » ee 
Rey, rato: TAR 
State | A?/A? (gs) State (a) A?/ A? (gs) 
(0, 1) 2-31 ons 2-50 (a/2)*-(1/2)40= 1-00 1-27 
(is 0) 3°95 eS 2-50 (SZ) (1i2)\= 0 0-40 
(2,0) | ~6 2 BD | 0-60 (3/2)? (1/2) | 0 0-04 
(eeb)e} eL 0% sed B) 0-40 (3/2)? .(1/2)?| 0 0-04 


§ 10. Discussion 


We have shown how the long f-lifetime of 14C can be reconciled with the 
spectrum of ™4N by using a tensor force which turns out to be small 
compared with the central and spin-orbit forces and in general agreement 
with that deduced from deuteron and triton data. Having done this in 
detail it is always instructive to look back in general terms at the solution. 
The wave function of the ground state emerged as 96% 1D, state, and 
the 2-31 mev level and ground state of “C are strong mixtures of *1S, 
and %8P, states. Hence the *C->14N f-decay can take place only theoneh 
the hall part (4%) of the ground state wave function, which is in the 
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3S, or 4P, states. In making the matrix cancellation of §3 we are 
setting conditions only on the small part of the wave function and we 
therefore need only a small tensor force. So far as the large part of the 
wave function is concerned the B-decay is L-forbidden, so that with no 
tensor force we should obtain log ft~5. The effect of this small tensor 
force on everything but the sensitive lifetime is therefore a perturbation 
on the usual intermediate coupling calculation and its inclusion in other 
nuclei will, in most cases, cause no sensible change in the intermediate 
coupling result. There may however be a few isolated cases, such as 
that considered here, where the tensor force is important. 

The inhibition of the M1 transition giving small A, which was discussed. 
in § 8, comes from a similar source since the large part of the wave function 
involves an S+D transition which is forbidden for M1 radiation. . 

Jancovici and Talmi (1953) pointed out that some confirmation that 
the large ft-value is produced by accidental cancellation is furnished by 
the much smaller value of log ft=7-3 reported by Sherr et al. (1955) for 
the 4O+'N decay which differs from the 4C+14N decay only in the 
Coulomb correction if we assume charge independence of nuclear forces. 
They showed that the Coulomb correction is just sufficient to destroy 
the cancellation, increasing the matrix element by a factor 10 and so 
reducing log ft from 9 to 7. A similar argument would apply to the 
calculations of this paper. 
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Note added in proof.—tlt has been brought to the author’s notice that 
Hird, Whitehead, Butler and Collie (1954) measured the branching ratio 
for the decay of the 3-95 Mev level in “N finding A=5+5%. This is in 
good agreement with the small value ~1°%, ‘ predicted ’ in § 8 of this paper. 


Further note added in proof.—Very recently Gove, Litherland, Almqvist 
and Bromley have also determined A through the reaction C!2(He®, py)N!4 
finding A=4-3-+-0-8°% in agreement with Hird et al. and with these calcu- 
lations. The fact that this value, although small, is about four times 
larger than the calculated value is consistent with the enhancement to be 
expected for an E.2 transition in this mass region from some weak surface 
coupling, see Elliott and Flowers (1955). 
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ABSTRACT 


It is a feature of many atomic diffusion jump mechanisms that successive 
jumps of an atom are ‘ correlated’ in the sense that the direction of a 
jump is not completely at random but is influenced to some extent by 
the direction of the preceding jump or jumps. General formulae are 
derived for the diffusion coefficient which take account of this correlation. 
The theory is applied to a calculation of the diffusion coefficient of solute 
atoms by the Johnson mechanism. 


$1. IyTRODUCTION 


THE starting point for a descriptive treatment of diffusion is the assump- 
tion of Fick’s law, which for one dimensional diffusion reads 


JanD) 15.29 . 27 


where J is the current density, dc/dxv is the concentration gradient and D 
is the diffusion coefficient. The purpose of this paper is to discuss an 
assumption which is frequently made in the derivation of the theoretical 
expression for D in terms of atomic quantities. We begin by summarizing 
the general steps in such a derivation to the point where this assumption is 
made. 

It is accepted that the fundamental process responsible for diffusion 
in crystalline solids is the continual migration of each of the atoms by 
an endless succession of jumps from site to site throughout the lattice. This 
entails a mixing of the atoms and so leads to the final effect of diffusion, 
viz., the removal of differences in composition. It is easily shown that a 
continual migration on the part of each atom does in fact produce a 
directed net flow down a concentration gradient at a rate proportional to 
it as required by the macroscopic eqn. (1.1). 

The usual and simplest way of doing this is to consider only those 
jumps that occur between two adjacent lattice planes (1 and 2) normal to 
the direction of the concentration gradient (see e.g. Cottrell 1955). 
J, which is simply the difference in the numbers of atoms jumping per 


+ Communicated by the Authors. 
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unit time in opposite directions between unit area of 1 and 2, is then 
found to be given by an expression of the form (1.1) with 


D=4I"s?, 


where J” is the probability per unit time that an atom will make any one 
of the possible jumps which take it out of the plane (1 or 2) and s is the 
interplanar spacing. The factor } arises because it is assumed that the 
jumps are random, in the sense that whenever an atom on one of the 
considered planes makes a jump the probability is } that it does so towards 
and on to the other plane and } that it jumps in the opposite direction. 

Now for a number of atomic jump mechanisms this is not true for every 
individual atom jump, and it is this assumption and the manner in which: 
the problem is to be handled when it is not valid that forms the subject of 
this paper. 

As a typical example consider diffusion by the vacancy mechanism. 
An atom which has just arrived on plane 1 from plane 2 by exchanging 
places with a vacancy, of necessity has this vacancy immediately adjacent 
to it on plane 2. This makes the probability greater than } that its next 
jump is back to 2, in which event the net contribution to J of these two 
jumps would be zero. If it arrived on plane | from the plane before 1, 
the probability is less than } that it contributes to J by jumping to 2, 
In either case the effect is to reduce J and hence D. To assume that the 
jump probabilities are always } would be to assume that immediately 
after a jump the vacancy responsible always migrated away from the 
atom and made further contact with it in only a random way, 
i.e., approaching from all directions with equal probability. In fact the 
probability of an atomic jump in one direction depends on the direction 
of the previous jump and the jump directions are said to be * correlated ’. 

' This correlation effect for vacancy self diffusion seems first to have been 

noticed by Herring (see Bardeen 1948, especially footnote p. 1404) and 
has been discussed by Bardeen and Herring (1952), who calculated for a 
quadratic layer lattice and for the body-centred cubic and face-centred 
cubic lattices, the decrease in D which it brings about. Since their treat- 
ment is rather condensed and since one of their equations (eqn. A.2) 
contains an elementary but important error we feel it worthwhile to discuss 
the topic again in simplified form. The corrected equations show that 
in both face-centred and body-centred cubic lattices D is reduced by 
about 20% 

There can be a much stronger correlation between successive Jumps in 
the case of the diffusion of solute atoms, at low concentration, by the 
Johnson (1939) mechanism. When there is an appreciable attraction 
between solute atoms and vacancies, comparatively stable solute atom- 
vacancy pairs will be formed and these can migrate as entities through 
the lattice. A particular jump of a vacancy can be with either its 
associated solute atom (probability per unit time, w,) or with any of the 
solvent atoms which are neighbours of both solute atom and vacancy 
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(probability per unit time for exchange with any one such solvent atom, 
w,). Ifw,>w, then there is a high probability that a second jump of 
a solute atom will be the reverse of its first jump. On the other hand if . 
w >We, then after one jump of a solute atom the vacancy will make on 
the average a large number of moves around the shell of nearest neighbour 
solvent atoms surrounding the solute atom, before exchanging again with 
the solute atom. By this time its position among the sites neighbouring 
the solute atom will be practically random. The degree of correlation 
can then be very large (wy >w,) or very small (w,>w,) depending on 
the relative magnitudes of w, and w,. We therefore expect D to contaim 
the rate of jumping not only of solute atoms but also of those of the solvent, 
so that the correlation here produces a qualitative change in the expression 
for D as distinct from the small quantitative correction that it makes in 
the vacancy self-diffusion case. Previous treatments of the Johnson 
mechanism are either incorrect or incomplete (Johnson 1939, Wyllie 
1947) and we shall therefore in §3 apply the general results of §2 to a 
calculation of D for this case. 

In the interstitialey mechanism of diffusion an atom in an interstitial 
position migrates by pushing a neighbouring atom into another interstitial 
site, itself occupying the normal site vacated. The path of a particular 
atom is then a sequence of jumps between alternately normal (nm) and 
interstitial (7) sites. A jump from an 7-site to an n-site is clearly uncor- 
related with previous jumps, the atom having an equal probability for 
jumping to any one of the allowed n-sites. However there is an appre- 
ciable probability that its next jump from the n-site will consist of being 
pushed back to its original 7-site by the interstitial atom created by its 
previous jump. In this example then only alternate pairs of jumps are 
correlated (i.e. pairs of the type i>n-+). 

There is clearly no correlation between successive jumps of a solute 
atom in an interstitial solid solution, for these jumps are from one inter- 
stitial site to another, each one completely at random. 

From these few but important examples we see that correlation is a 
very frequent but not a universal feature of atomic jump processes in 
crystals. 


§2. THE GENERAL FORMULA FOR THE DIFFUSION COEFFICIENT 


When correlation effects exist it is clearly inadequate to consider only 
those jumps that occur between two adjacent planes of atoms for we 
require to know the directions from which atoms arrive on these planes, 
and are therefore interested in jumps made from other planes on either 
side of them. This suggests that we could take account of correlation if 
we considered jumps between larger numbers of consecutive planes. One 
of us has in fact recently treated the Johnson mechanism by this method 
and the results agree with those of §3 (Lidiard 1955, particularly § 4). 

A more satisfactory method, yielding more readily to general treatment 
is to consider not merely jumps between a given number of planes but 
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the whole sequences of jumps that constitute the paths of atoms migrating 
through the lattice. Such an approach to diffusion problems is known as 
the ‘random walk method’. In the present case it is more correct to 
speak of correlated random walks. 

If successive jumps are of vector magnitude r,, r,, etc. then after 
a time ¢ a total of n=TJ¢t will have been made + and the atom will be a 


n 
vector distance R(é)= 2 r, from its original position. R2(t), the average 
i=1 
of R(t) over a large number of identical atoms is simply related to D; 
it can be shown that for one-dimensional diffusion 


dc R*(t) dc 


ot BECO? ea) 
(see e.g. Joos 1934, §7 of chapter XX XIV). 
Whence, 
ee) 
D=— (2.2) 


The next step in the argument is to show that R%(t) is proportional to 
¢ and to relate it to the individual jumps which compose R. We have 


Ser ei rer ya ta a Sf. * (2:3) 


We shall assume for simplicity that the individual jump vectors are of 
equal magnitude vr. This is true for the commonly occurring face-centred 
and body-centred cubic lattices and also for the hexagonal close packed 
lattice. Since the average value of r;.r,,; will be independent of 7 in 
the case of vacancy diffusion we can write (2.3) in the form 

—— nm—1 BS a ier, 

R2(é)=nr2+ 212 2X (n—j) cos 0;(vacancy). . . . . (2.4) 

j=1 
where cos 0; is the average value of the cosine of the angle between the 
ith and the (¢-+7) jump of an atom. rhs 
If the direction of one jump is completely independent of the direction of 


the preceding or earlier jumps the jumps are uncorrelated and cos 6;=0 


for all j, since all possible values of 0; are equally probable. R*(Z) is then 
proportional to ¢, since n=J%, and D becomes 


ae eee eee ee. (2.5) 


This is identical with the previous expression for ) when the relations 
between I’ and I” and between r and s are inserted. 
; 
+L is the total number of jumps which an atom makes per unit time and. is 
not to be confused with Ir’ introduced earlier. The relation between the two 
quantities is easily determined once the lattice is specified. 
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Jan 


As we have noted in $1 (2.5) is applicable to interstitial diffusion in 
interstitial solid solutions, but is not exactly applicable to vacancy or to 
interstitialey diffusion. However it is to be emphasized that so long as 
we are dealing with self diffusion in a pure crystal eqn. (2.5) is, 
as we shall see, in error only by a numerical factor not very different from 
unity, and is not qualitatively incorrect ; most of the theoretical analyses 
which have been made of D therefore still stand (see for example Zener 
1952). In eases like that of the Johnson mechanism, (2.5) is considerably 
in error and the correct equation differs from it radically. 

To discuss correlation effects we return to eqn. (2.4) and consider the 
quantities cos@;. Firstly in order to see the relation between cos 6; 
and cos 6, , consider the vectors r,,r;,;-; and r,,; Their poles in 
stereographic projection form the corners of a spherical triangle the sides 
of which are the angles 6,, 6;_, and 6;. Writing down the standard 
relation} between these angles, and taking averages on both sides we get 


cos 0; = cos 4; ,cos@,-+ sin@,,sin@,;cosé . . . . «+ (2.6) 


where ¢ is the angle between the two planes defined by r; and r,,,_, and 
by r,,;andr,,, The averages of the product terms on the right-hand 
side can be replaced by the products of the averages of their separate 
factors, provided these factors are statistically independent of one another. 
This condition is satisfied for vacancy diffusion since cos 6,, the mean 
cosine of the angle between two successive jumps, is independent of the 
jumps which have gone before. In other words there is a direct correla- 
tion only between consecutive jumps, the correlation between one jump 
and a jump later than the next, being indirect and determined only by the 
direct correlations between successive intermediate jumps. Also the 
probability of a given angle 6 between two jumps depends only on @ and 
not on the orientation of the plane of the jumps—at least for lattices where 
all atom jumps are of equal length and of equal probability, so that 


cos¢=0. Thus (2.6) yieldst 


cos. 0, = cos 0,.,..cos 6, = (cos @,)), ... 8 2 {2a 


Substituting (2.7) into (2.4) and passing to the limit n, t+ oo we get 


D={e ED, (vacancy) Ak ei hive eng cn meee 


1 See e.g. Introduction to Crystallography by F.C. Phillips, p. 182. 

~ A more rigorous derivation of eqn. (2.7) is possible. It has been shown 
(RK. Howard and A. B. Lidiard) that (2.7) is valid for vacancy diffusion in all 
lattices having one atom per unit cell and all jump vectors of the same magni- 
tude. It is also true for the plane honeycomb layer and the diamond lattice. 
The apparent contradiction of (2.7) provided by Bardeen and Herring’s eqns. 
(A.3) vere (A.5) for the plane quadratic layer is not real since their eqn. (A.5) is 
incorrect. 
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where we have written C for cos cos 6, the average of the cosine of the angle 
between consecutive jumps. The term (1+-C)/(1 
as the ‘ correlation factor ’. 

By a detailed calculation (the principles of which we illustrate in §3 when 
we consider impurity diffusion by the Johnson mechanism) Bardeen 
and Herring obtained numerical values for C for self diffusion in body- 
centred and face-centred cubic lattices. Their results are respectively 
O(b-c.c.)=—0-124 and C(f-c.c.)=—0-109. Substituted into (2.8) 
these values give correlation factors of 0-78 and 0-80 respectively. 

The derivation of (2.8) assumes that the direct correlation between a 
pair of consecutive jumps is the same for any such pair. We have already 
seen that this is not true for interstitialcy diffusion where only alternate 
pairs of jumps are correlated at all. This simplifies considerably the 
derivation of the correlation factor. Of the terms r,.r,,; in eqn. (2.3) 
all those with 7>1 are zero, and of those for which j=1 only alternate 
ones are non-zero. For this case therefore 


R2(t)=nr2 (1--cos 0 0) 


—C') we may refer to 


so that 
D=}{Ir?(1+ cos @) (interstitialcy). » Eee eee ke) 


Here cos @ is the mean value of the cosine of the angle between consecutive 
in and ni jumps. 


§3. Lwpurtry DIFFUSION BY THE JOHNSON MECHANISM 


The Johnson mechanism of the diffusion of a solute in a dilute alloy 
assumes the existence of short range forces of attraction between a solute 
atom and a vacancy in the host lattice. In the following discussion we 
shall assume that we are dealing with a face-centred cubic lattice and that 
the attraction between an impurity atom and a vacancy may be neglected 
unless they are nearest neighbours. The average fraction of solute atoms 
with which a vacancy is associated will be denoted by p. The number of 
times that the associated vacancy and impurity change places with each 
other per unit time is w,; w, is the frequency with which an associated 
vacancy will jump from one position to another particular nearest neigh- 
bour position of the impurity. (Each pair of nearest neighbours in a face- 
centred cubic lattice has four common neighbours). Hence ¢J/7? is 
seen to be equal to 4 a2w,p when the distance between nearest neighbours, 
r, is put equal to 4/2a. 

It remains to calculate C. Let us do this first for the limiting case where 
the vacancy-solute pairs dissociate very much less frequently than they 
change their orientation. Suppose that the initial jump of the impurity 
was from position} (011) to (000) so that we start with the impurity at 
(000) and the vacancy at (011). We wish to calculate the average cosine 
of the angle between the vector (0—1—1) and the next impurity jump 


— 


+ For convenience we omit the factor a. 
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vector. The first jump of the vacancy at (011) is back to (000) with 
probability w,/(4w,-++w,). This means that the impurity jump vector 
is (011), ie. opposite to (0—1—1) so that the cosine of the angle is —1. 
The contribution to the average cosine C from this possibility is thus 
(—1)w,/(4w,-++-w,). On the other hand it is easily verified that the first 
vacancy jump could have been to one of the positions (110), (101), (—110) 
and (—101) each with probability w,/(4w,+-w.). We now consider the 
second vacancy jump. First we see that we can forget those systems in 
the ensemble where the vacancy went to (000) since we are interested only 
in calculating the mean cosine of the angle between two consecutive jumps 
of the impurity. Thus we consider that the second jump of the vacancy 
may be from any of the above four positions to (000), each with probability 
W,/(4w,-+w,). The cosine of the angle is —4} for each of the four positions. 
There is therefore a contribution to C from the second vacancy jump of 
amount, 
(act) We 4wy, ; 

“" (4w,+W,) (401+ Wg) 
By following the vacancy around the twelve neighbouring sites to (000) 
in this way, we could write down a contribution to C coming from every 
vacancy jump. By induction we can then write down and prove the 
nature of the term in C' coming from the n™ jump. Summation of the 
series for C so obtained gives 


C= =Wyf(2 05) SS ee 


However a more elegant formulation is possible and is necessary for 
more complicated cases than the one we are considering in this section. 
After the nt vacancy jump let the probability of occupation of the site 
s by the vacancy be p,(s). Thus with the above initiai conditions 
p(011)=1 and all other p(s) are zero, whilst p,(110), p,(101), »,(—110) 
and p,(—101) each equal w,(4w,+wy,). By ordering the sites s in some 
suitable manner the various p,,(s) coefficients can be written as a column 
matrix, p,. The contribution of the n+1‘ jump to C can then be 
expressed as the matrix product 7p, where 7 is arow matrix each of whose 
elements is the probability, w,/(4w,+-w,.), that a vacancy will jump to 
(000) multiplied by the cosine of the angle for the appropriate site. If 
Pp, has been contracted because some elements are always equal to one 
another by symmetry, then 7 must also be contracted and each element in 
addition to the above two factors will contain a third factor which is the 
number of distinct sites which are represented by the single probability 
corresponding to that element of r. We have then 


Casa Tia ee Ce oe a een) 


Also the probability coefficients p,4,(s) are related to the p,(s) by a 
set of linear equations. So we write 


Pnsi=Ppy | 
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where P is the matrix of the coefficient in these equations. Hence 
0 
Cred Pp): ee eee os 8 (3,3) 
n=0 
A glance at fig. 1 shows that in the present case, owing to the symmetry 
Rigel 


X = IMPURITY ION 


C= NEAREST NEIGHBOUR 


The twelve nearest neighbours to an impurity atom (x) in the face-centred 
cubic lattice. 


around the axis (011), there are only five distinct probability coefficients. 
These we take to be (1) p,(011), occurring once, (2) p,(110), occurring 
four times, (3) p,(01—1), occurring twice, (4) p,(10—1) occurring 
four times and (5) p,(0—1—1) occurring once. With these coefficients 
arranged in this order in a column matrix we easily find that 7 must be 
the row matrix 
| | mee 02 ion es (8.4) 
where we have used g as an abbreviation for w./(4w,+w,). By writing 
down the linear equations connecting p,,,, with p,, we easily find that P is 


Om 4s 30.0 


Th hk a 
P=f |0 2 0 2 0 (3.5) 
Op aa ie elar |: 
Om 0 0254=— 0 
where f=w,/(4w,+w,). From (3.4) and (3.5) we find 
) 7P=2fr, 


whence 
Tha (2 er. 
Finally, since py is the column matrix 
po={l 0 0 0 Oh, . 
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we see that (3.3) gives 
C=7p,. 2 (2f)"=—g-2f)"*; J a ee oe 
n=0 


which by insertion of the expressions for f and g leads to 


O= —weal(Qwy+w,), 2 eee 
as stated above. 

Up to this point we have assumed that dissociation jumps by the 
vacancies attached to impurity ions are very rare. It is a simple matter 
to remove this restriction if we assume that the return probability of a 
dissociating vacancy may be neglected in calculating C. Let 4, be the 
probability per unit time that an associated vacancy will make one 
particular dissociation jump, i.e. a jump to a particular site which is not 
a nearest neighbour of the impurity. An associated vacancy has seven 
such jumps open to it. The preceding calculation remains unaltered 
except that the relative probabilities f and g are now respectively 
W,/(Wot4w,+7k,) and w,/(w.+4w,+7k,). Use of eqn. (3.6) then gives 

‘C= — Wo] (2, twet Tee 2 aa Sa ee 

If we pass to the limit when w,=w,=h,, i.e. when the ‘impurity ” 
is a tracer atom, then (3-7) predicts that C=—0-10. This may be com- 
pared with the Bardeen and Herring figure of 0-109 for a face-centred 
cubic lattice. The difference arises from the fact that we have 
neglected the return probability of a ‘ dissociating’ vacancy. Bardeen 
and Herring included this probability by taking explicit account of the 
probability coefficients p,(s) for next nearest neighbour and more distant 
sites. The principles of the calculation are the same although the detailed 
computation becomes rapidly more involved as the order of P increases 
especially when 7 is not an eigen vector of P as it is here. It is however 
certain that the error in (3.7) decreases as the binding between the impurity 
atom and the vacancy becomes ROR 2ST, i.e. as k, decreases relative to 
w, and wy. 

Insertion of (3.7) into (2.8) gives the final expression for the diffusion 
coefficient appropriate to the Johnson model in a face- centred cubic 
lattice, 

p—wP wW,(w,+ Tk, /2) 
3 (wy +W.+ 7k, /2) 
—in agreement with the result which was previously obtained using kinetic 
theory (Lidiard 1955). 


(3.8) 


§4. SUMMARY AND CONCLUSION 


In the present paper we have given a simplified account of the random 
walk method of dealing with correlation effects in atomic diffusion in 
crystals. The general formula (2.8) for the coefficient of diffusion by 
the vacancy mechanism in a lattice where the elementary jumps are of 
equal length, expresses the fact that direct correlations exist between 
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consecutive atom jumps and that they exist equally between every pair 
of consecutive jumps. When the institialey mechanism is operative 
such correlations exist only between pairs of type in followed by ni, 
and the diffusion coefficient is then given by (2.9). In §3 we illustrated 
the matrix formulation of the method of evaluating C by a calculation of 
the Johnson model of impurity diffusion in a face-centred cubic lattice. 
The results agreed with those obtained previously using a kinetic method. 
At present we are investigating by a random walk method the self diffusion 
of solvent atoms by vacancies which are strongly bound to impurity 
atoms. A practical example of this state of affairs is the diffusion of 
radioactive sodium in a erystal of sodium chloride containing divalent 
impurity cations, e.g. CdCl, (Etzel and Maurer 1950, Aschner 1954). 
A reliable kinetic treatment of this problem would be extremely 
cumbersome. 

Note added in proof. The authors are indebted to Prof. C. A. Coulson 


for pointing out that equation (2.8), or rather the equation for R*(t) 
corresponding to it 


R?2(t)=nr?(1+C)/(1—C), 


has been employed in the theory of rubber elasticity for a number of 
years. It gives the distribution of lengths of a chain of n carbon-carbon 
bonds when there is free rotation around each bond but the bond angle 
is constant (6,). See e.g. L. R. G. Treloar, The Physics of Rubber 
Elasticity (Oxford: Clarendon Press), 1949. 
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ABSTRACT 


A differential method has been used to measure the changes in density 
which accompany the annealing of plastically deformed nickel. The 
results of these measurements are correlated with earlier measurements of 
stored energy, electrical resistivity and hardness for the same material. 
There are two sudden increments in density ; the first is attributed to the 
disappearance of vacancies and the second to recrystallization. These 
are superimposed on a gradual increase in density attributed to recovery. 
There is a considerable discrepancy between the densities of dislocations 
calculated from the stored energy measurements and from the density 
measurements. 


§ 1. INTRODUCTION 


In a recent paper (Clarebrough, Hargreaves and West 1955, hereinafter 
called I) the release of stored energy during the annealing of plastically 
deformed nickel was correlated with changes in hardness and electrical 
resistivity. The conclusion was reached that the disappearance of both 
vacancies and dislocations contributed to the release of energy. In an 
attempt to separate the contributions of these two types of defect, 
measurements have been made of the changes in macroscopic density 
which take place during annealing of plastically deformed nickel. 

Preliminary results, given in I, were obtained by Mr. G. A. Bell of 
the National Standards Laboratory, C.S.I.R.O., Sydney. He measured 
the absolute density of a specimen as deformed, again after heating to a 
temperature considered high enough to allow the vacancies to disappear 
and again after further heating to a temperature in excess of the recrystalli- 
zation temperature. The details of the method used will be published 
elsewhere (Bell 1956). These preliminary results did not give sufficient 
detail of the manner in which the change of density was related to the 
temperature of annealing, and more detailed experiments were required. 
These are described in the present paper. A differential method for 
measuring small changes in density was developed simultaneously by 
Mr. Bell and by the authors. This method has been employed in making 
the measurements described here. 


+ Communicated by W. Boas. 
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§ 2. EXPERIMENTAL 
(a) Preparation of Specimens 


The nickel used for the density measurements was of commercial purity 
(99 6°, Ni) and from the same batch (K) as that used for the measure- 
ments of stored energy, electrical resistivity and hardness (I). Deforma- 
tion in torsion and in compression was carried out on specimens of the 
same dimensions and in the same manner as described in I. The bars 
deformed in torsion were ground lightly to remove the roughened surface 
produced by the heavy deformation. Specimens for the measurement 
of density, 13m. in length and ? in. in diameter, were then cut from the 
ground bars and the end faces ground flat. Since the bar deformed 
by compression was not in a suitable form for grinding, the specimen for 
the measurement of density was produced by turning. All cutting, 
grinding and turning operations were done under a copious flow of coolant. 

One of the deformed specimens was annealed by heating at a rate of 
6°c per min to 900°C and this specimen was used throughout the experi- 
ments as the standard to which all deformed specimens were referred. 
In each case the weight of the deformed specimen in air was made equal 
to the weight of the standard in air within | mg. 


(b) Procedure 


The differential method involved determination of the difference in 
upthrust between the standard specimen and a deformed specimen 
when both were immersed in a liquid of high density. 

An analytical balance with a capacity of 200 g and a nominal sensitivity 
of 1/10 mg was used in the experiments. A concave mirror was attached 
to the beam of the balance to focus a beam of light onto a scale at a 
distance of 5 metres ; rest points were determined on this scale. 

The balance was mounted rigidly above a water-bath and the tempera- 
ture of this bath was controlled at 27-0-01°c. A glass vessel containing 
the working liquid (1, 1, 2, 2 tetrabromo-ethane) was supported centrally 
in the water-bath. This vessel is described elsewhere (Bell 1956) and the 
design allows the liquid to be stirred without introducing a mechanical 
stirrer. The stirring was stopped during weighing. For weighing in 

iquid the specimens were held in nickel wire cages, which in turn were 
suspended by fine nickel wires from hooks underneath the balance pans. 

Before making a measurement the cages were suspended in the working 
liquid and sufficient time allowed for the cages and the liquid to reach the 
operating temperature. The rest point of the balance was then determined 
with the standard specimen and the deformed specimen on the pans and 
the cages in the liquid. The specimens were then transferred to the 
cages in the liquid and, after they had reached the temperature of the 
liquid, the rest point of the balance was determined again and adjust- 
ments made with a rider to restore the rest point to approximately that 
for the specimens in air. The specimens were then removed from the 
liquid, washed in alcohol, dried and the rest point in air was checked. In 
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this way the difference in upthrust between the deformed and the 
standard specimen was measured. 

The difference in upthrust was measured for each specimen after 
deformation and after heating to various temperatures at a rate of 6°C 
per min. From the differences in upthrust measured before and after 
any heating the corresponding fractional change in density may be 
calculated knowing the density of the liquid. This was determined from 
time to time by weighing in air and in the liquid a standard nickel specimen 
whose absolute density was known accurately from Bell’s measurements. 

For heating, the specimens were placed in a nickel boat inside a silica 
tube which was evacuated and then heated in a furnace to a particular 
temperature, withdrawn and allowed to cool. The vacuum was then 
broken, the specimen removed and the measurement of density made. 
The specimen was then replaced in the tube which was again evacuated 
and then returned to the furnace held at the particular temperature and 
heated to the next temperature required. 

The differential method has several advantages over the measurement 
of absolute densities for determining small changes in density. The change 
in density is measured directly instead of being obtained as a small 
difference between two large measured quantities. Since two specimens 
of very nearly the same volume are used any changes in the density of 
the working liquid only cause second order errors in the result. For this 
reason the temperature of the working liquid need not be controlled 
within the very close limits required for the measurement of absolute 
density. The results of the differential method are not sensitive to 
errors in the balance nor do complications arise from the precision of 
weights or the use of buoyancy corrections. 


§ 3. RESULTS 


Figure 1 shows the fractional change in density 4D/D and the power 
difference 4P as functions of the temperature for specimens deformed in 
torsion (nd/l= 2-34, where n=number of turns, d=diameter and /=: gauge 
length) and heated at a rate of 6°c per min.: The power difference 4P 
represents the rate of release of stored energy (Clarebrough, Hargreaves, 
Michell and West 1952). The density results were obtained with five 
separate specimens and the 4P curve is the mean of four separate 
determinations. 

The influence of the degree of deformation on the changes in density is 
shown in fig. 2 where 4D/D is plotted as a function of temperature for three 
deformations in torsion corresponding to ndil values of 0:47, 1-41 and 2-34. 

Figure 3 shows 4D/D as a function of temperature for a specimen 
deformed approximately 80° in compression. It can be seen that this 
curve has the same form as the 4D/D curves in fig. 2 for specimens 
deformed in torsion. Thus the type of result is not influenced to any 
extent by the inhomogeneous deformation in torsion, as was also found 


for the determinations of stored energy, so that the results should be 
generally applicable. 
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Fig, 1 
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The rate of release of stored energy 4P and the fractional change in density, 
measured at 27°c, 4D/D, as functions of temperature for nickel deformed 


in torsion to nd/]=2-34. 


Fig. 2 


60 


52) 


40 


20 


Meen of 5 
specimens 


° exe) 200 300 400 500 600 700 800 900 


TEMPERATURE °C 


The fractional change in density, measured at 27°C, as a function of tem 
perature for nickel deformed to various extents in torsion. 
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The fractional change in density, measured at 27°c, as a function of tem- 
perature for nickel deformed 80° in compression. 
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The total change in density, curve A, and the portion attributed to vacancies 
curve B, as functions of the amount of deformation in torsion. 
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For all deformations the 4D/D curve indicates a gradual rise in density 
on which two sudden increments in density (AB and CD in fig. 2) are 
superimposed. In fig. 4 the total fractional change in density together 
with the fractional change involved in the first increment, AB, are plotted 
against the torsional strain. The values for the first increment were 
obtained by extrapolating the portions BC of the curves back to 200°c. 
This procedure can only give a very approximate result for reasons 
discussed below, but the trend of the curve is clear. 
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A typical result from (I). The rate of release of stored energy 4P, the incre- 
ment in electrical resistivity measured at 0°c, dp, and the hardness 
measured at room temperature, V.H.N. as functions of the temperature 
for nickel deformed in torsion to nd/l=1-87. 


§ 4. Discussion 


The changes in density observed should be considered in relation to 
the determinations of stored energy, hardness and electrical resistivity 
which had been made previously on the identical material. A typical 
set of these results from I is shown in fig. 5. The release of energy corre- 
sponding to the first peak in the 4P curve, area a’, is accompanied by a 
decrement in the resistivity without any change in hardness and therefore 
has been attributed to the disappearance of vacancies produced by 
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deformation. On a vacancy inodel Nicholas (1955) obtained good 
agreement between calculated and experimental 4P curves. 

Now it is clear from fig. 1 that the first sudden increase in density 
AB takes place in the same range of temperature as the release of energy 
corresponding to area a’, the first peak in the 4P curve. It was shown in 
I that this range of temperature was independent of the amount of 
deformation and the curves in fig. 2 show that this is also true for the 
range of temperature over which the first sudden increase in density 
occurs. Further, the amount of energy represented by area a’ (fig. 13 
in I) and the magnitude of the first increment in density (curve B, fig. 4) 
are similarly almost independent of the amount of deformation. Thus 
the conclusion reached in I that vacancies created during deformation 
disappear during this range of temperature is supported by the present 
observations of the corresponding increase in density. 

In estimating the energy of a vacancy it was assumed previously that 
the disappearance of vacancies was responsible for the whole of the increase 
in density observed, in the preliminary measurements, up to the tempera- 
ture at which the release of energy represented by area a’ was complete. 
This assumption is not supported by the present results, as there appears 
to be a general upward trend in the 4D/D curve on which the sudden 
increment is superimposed, and the extrapolation referred to in §3 must 
be carried out to find the change im density associated with area a’. 
The proportion of vacancies is therefore less than that calculated on this 
assumption and it now appears that the estimate of the energy of vacancy 
as 1-6 ev must be increased. However, the increment in density due to 
the disappearance of vacancies cannot be determined exactly from the 
present results for the following reasons. 

Firstly, the necessary extrapolation cannot be made accurately without 
measuring the increase in density at smaller intervals of temperature. 
Secondly, the technique of annealing in vacuo adopted in this work 
requires an appreciable time for a specimen to reach the end point of the 
previous annealing, before heating through the next range of tempera- 
ture, and also the rate of cooling after each heating is quite slow. A 
specimen on which the change in density is determined is therefore at 
elevated temperatures for a much longer period than a specimen used for 
a determination of the 4P curve by uninterrupted heating in the calori- 
meter. As it has been found that the temperatures at which the peaks in 
the 4P curve occur are somewhat dependent on the rate of heating, an 
exact correspondence of the 4P and 4D/D curves, with respect to tempera- 
ture, is not to be expected. Further, the closer the intervals of tempera- 
ture at which the increase in density is determined the more serious is 
this discrepancy likely to become. 

It is hoped to resolve this difficulty by improving the technique of 
annealing and making experiments on a much purer grade of nickel. 
Preliminary experiments with this material indicate that the increase 
in density due to the disappearance of vacancies may be more clearly 
resolved, 
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The release of energy corresponding to area (a”+B) in fig. 5 is accom- 
panied by only a slight decrease in hardness and a gradual decrease in 
electrical resistivity, on which is superimposed the more sudden decrease 
due to vacancies. This is very similar to the gradual increase in density, 
which takes place in the corresponding range of temperature, on which the 
increment due to vacancies is superimposed. The release of energy 
accompanying these changes has been attributed in I to recovery involving 
the annihilation and re-arrangement of dislocations. If it is accepted that 
the presence of dislocations in a metal reduces its density this gradual 
increase in density may be attributed to the operation of these same 
processes. 

The second sudden increase in density, CD, takes place in approximately 
the range of temperature in which the second peak in the AJP curve, 
area C, is observed. Measurements of hardness and resistivity such as 
those illustrated in fig. 5 together with metallographic and x-ray observa- 
tions show that this peak corresponds to recrystallization (I). Thus, 
the second sudden increment in density may consistently be attributed 
to the decrease in the density of dislocations during recrystallization, 
assuming again that dislocations cause a reduction in density. The 
AD/D curve begins to rise at a temperature rather lower than that at 
which the peak, area C, in the 4P curve begins. This is probably due to 
the deficiencies in the technique of annealing, which were pointed out 
above, since the total time durmg which the specimen has been at 
elevated temperatures has become considerable at this point in the 
experiment. Despite this effect it is clear from fig. 2 that the sudden 
increase in density, CD, attributed to recrystallization, occurs at higher 
temperatures for smaller amounts of deformation in the same way as do 
the relevant peak in the 4P curve and the accompanying changes in 
properties. 

It is of interest to compare the density of dislocations obtained from the 
results of the measurements of stored energy given in (I) and from the 
results of the present measurements of density. The appropriate values 
for two degrees of deformation are given in the table. 
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In calculating the density of dislocations from the energy measurements 
it is assumed that the liberation of energy in area (a’4+B-+C) is due to 
a decrease in the density of dislocations and that the energy of a dislocation 
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in nickel is 7x 10-4 erg em=!. As pointed out in (I) the first assumption 
is not fully justified as some interaction energy is almost certainly involved 
in the liberation of energy in area (a’+B-+C). In addition, the value 
assumed for the energy of a dislocation is somewhat arbitrary and a 
value as low as 2-5 10-4 erg em~+ might be taken for dislocations in 
heavily cold worked nickel. The corresponding densities of dislocations 
would then be 5:9 x 10! and 8-8 x 1044 lines cm. 

In calculating the density of dislocations from the results of the present 
measurements of density, the total change in density, less the contribution 
due to vacancies, is assumed to be due to a decrease in the number of 
dislocations. It is assumed that an edge dislocation causes the same 
change in density as a row of vacancies of the same length, which is 
probably a reasonable upper limit, and that half of the dislocations are 
of the screw type, which make no contribution to the change in density. 
On this model, if the dislocation lines run between nearest neighbours, 
distant a/4/2, the total density of dislocations V which causes a fractional 
change in density 6 is N=4y/26/a? lmes cm”. 

Taking the lower limit for the energy of a dislocation and this upper 
limit for the density change due to a dislocation there is still a considerable 
discrepancy between the densities of dislocations calculated from the 
energy measurements and from the density measurements. Thus it 
does not appear possible to explain the present results on the basis of the 
model suggested, involving only dislocations, without invoking still lower 
values for the energy of a dislocation and consequently densities of dis- 
locations considerably greater than 10? lines cm~?. Such high densities 
of dislocations do not appear probable, on the available evidence, and 
thus it appears, that either the change in density due to a dislocation is 
higher than we have supposed, or that some other type of defect con- 
tributes to the change in density. 

A similar anomaly has been found by Van Bueren (1955) for the increase 
in electrical resistivity due to deformation in copper and he has suggested 
that point defects are responsible for most of this increase, though it is 
not clear by what mechanism these are stabilized at relatively high 
temperatures. ‘ 
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ABSTRACT 


A series of experiments, based on electrical resistivity measurements, 
have resulted in the following conclusions: The activation energy for 
formation of vacancies is about 1-0 ev and for migration about 0-9 ev. 
The ordering process, involving the interaction of vacancies with wrong 
atoms, is complex ; at least two relaxation processes are shown to operate. 
The decay time of an excess vacancy concentration appears to be limited 
by the presence of vacancy capturing impurities, the concentration of 
which can be varied by suitable heat treatment. Interstitials are more 
mobile than vacancies at a given temperature, but these too are probably 
trapped and can in turn contribute to the limitation of the vacancy decay 
time. 


$1. INTRODUCTION 


THE present work on the alloy Cu,Au originated as a study of the dis- 
ordering effect of neutron bombardment. However, when an ordering 
effect was also discovered interest turned in this direction (for a review of 
these effects see Kinchin and Pease 1955). Ordering requires the micro- 
diffusion of copper and gold atoms which, as with diffusion in face-centred 
cubic metals generally, is believed to occur through the migration of 
vacancies. Thus the objective has become, primarily, the study of the 
behaviour of vacancies in the alloy. Although it has not been possible to 
substantiate the vacancy mechanism directly the self consistency of the 
hypothesis is well supported. 
The general technique was to study changes in the electrical resistivity 
(always measured at 0°c) brought about by various treatments including 
bombardments and quenches. By this technique the alloy appears to be 
so sensitive to small concentrations of vacancies that bombardments 
with y-rays (discussed in the appendix) and low temperature quenches 
(from ~300°C) were convenient. | 

The description of the experimental work begins with preliminary 
experiments in support of the vacancy hypotheseis. Some plausible 
assumptions concerning the vacancy mechanism are then made and these 
are applied to determine H, the activation energy of formation, to study 
the interaction with wrong atoms and to find the factors which control 
the decay time of excess vacancies. In addition, some attention is given 


to the behaviour of interstitials. 
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§ 2. EXPERIMENTAL TECHNIQUE 


Much of the experimental technique has already been described 
(Dugdale and Green 1954—henceforth referred to as P1) and so only 
those aspects of particular relevance to the present work are mentioned 
here. 

The alloy was made by Johnston, Matthey and Co., Ltd., in the form 
of wire of diameter 0-025 mm, from high purity copper and gold. An 
analysis, made by the Chemistry Division at Harwell, gave 25:8--0-2 
atomic °%, gold, with silver, iron and silicon as major impurities with 
concentrations between 0-01 and 0:1 atomic °% 

Specimens were made from portions of the wire subjected to an initial 
heat treatment of 2 hours at 900°c, in an argon atmosphere, to remove 
coring. After suitably mounting on mica, they were given a further 
treatment at 600°C to remove any cold work effects. In the early stages 
a half hour treatment was given, but this was later increased to four hours. 
All the specimens mentioned in § 4 et seg. received the longer treatment. 

The majority of the experiments were concerned with the state of long 
range order. Normally, specimens were prepared in this state by heating 
for periods of 5 to 8 days at 370°C, after which they were cooled to room 
temperature in a period of about 5 minutes giving a value of the electrical 
resistivity p, at 0°c, of about 642 cm (aocording to Siegel 1951 domain 
growth occurs most rapidly at ~370°c). It was found that after such 
treatment (henceforth a specimen prepared in this manner will be termed 
‘partially ordered ’) equilibrium at lower temperatures (>250°c) could 
be closely approached in relatively short times; fig. 2.1 shows, for example, 
the relaxation towards equilibrium at 250°c and 290°c. Conversely, for 
‘specimens used in experiments during which the degree. of order was 
increased, the partially ordered state could be recovereed by a short heat 
treatment at 370°C of only 2 hours. 

Values of the equilibrium resistivity pe have been obtained, by quenching, 
from several specimens during the experiments; these are plotted in 
fig. 2.2 as a function of temperature. Also shown are the values for a 
domain size of ~10~° cm deduced from the work of Sykes and Jones (1938) 
who used an alloy very close to the stoichiemetric composition. It can 
be seen that the values obtained from the present alloy are very reasonabie. 
“Co y sources, giving uniform fluxes at the specimen of from 1014 to 1015 
y em? hour~!, were used for irradiations. 


§ 3. PRELIMINARY EXPERIMENTS 


The quenching in of vacancies is conveniently demonstrated on a 
tempering curve. For this purpose, a partially ordered specimen was 
quenched from 380°C and then heated at progressively increasing tem- 
peratures for 1 hour periods. The resistivity-temperature curve so 
obtained is shown in fig. 3.1 (a). The general shape of the curve, i.e., the 
passage through a minimum at ~300°c, is characteristic of the alloy, but 
the important aspect, for the present purpose, is the inflection in the 
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neighbourhood of 170°c. The rate of ordering is noticeably greater just 
below this temperature than it is just above ; in fact, the ordering rate at 
160°C is not again equalled until the temperature exceeds 200°C. This 
behaviour is just what might be expected if excess vacancies, frozen in 
by the quench, are mobile and decay to equilibrium at ~160°c. The 
higher ordering rates above 200°c are attributed to the thermal formation 
of sufficient vacancy concentrations. 

Figure 3.1 (6) shows a similar tempering curve taken after the specimen 
was cooled from 380°c to room temperature relatively slowly (cooling 
time ~5 min) and then heated at 160°c for 16 hours. No ordering is 
noticeable until a temperature of ~200°c is reached, which is consistent 
with the above interpretation. 

Quenching from lower temperatures gives smaller but similar in- 
flexions in the neighbourhood of 170°c. Such quenches, of course, also 
result in a lower degree of order and initial resistivity. Quenching from 
425°c (37°c above the critical temperature), on the other hand, results in 
the disordered state with a high resistivity. A tempering curve for this 
case (fig. 3.1 c) shows hardly any change in resistivity until ~250°c is 
reached. Although it has been shown by Skyes and Jones (1936) that 
considerable ordering does occur after such a quench it is evident that the 
resistivity is not sensitive to it. A state of long range order is therefore 
essential for the sensitive detection of excess vacancies at low temperatures. 

Low temperature ordering may be induced in a slowly cooled specimen 
by means of a bombardment with Coy rays. The tempering curve (a) 
in fig. 3.2 shows the behaviour before bombardment (this time with half 
hour periods at each temperature), and curve (b) after bombardment with 
1-5 x 1016 y-rays cm~?. 

The isothermal decay of excess vacancies at 160°c and 190°c, after 
bombardement with 1-5 x 10!® y-rays cm~?, is shown in fig. 3.3. The first 
25 hours of curve (a) shows the small but detectable thermal ordering at 
160°c. The bombardment then produced an acceleration in ordering rate 
which decayed in about 3 hours. At 190°c (curve (b)) the thermal ordering 
rate is appreciably higher and the decay time of the excess vacancies is 
reduced to about | hour. 

An attempt was made to observe the converse of this process at 190°c, 
i.e. the growth to equilibrium of an initially too small vacancy concentra- 
tion. The previous experiment, involving an increase of temperature 
from 160°c to 190°c and hence an increase in the equilibrium vacancy 
concentration, suggested that it would be a very small effect. A partially 
ordered specimen was heated at 190°c for two hours to reach the equi- 
librium vacancy concentration. The ordering rate as a function of time 
during a subsequent hour is plotted in the first part of fig. 3.4. The 
specimen was then heated for long periods at 130°c and 100°c (respectively 
24 and 134 hours) to reduce the vacancy concentration as much as possible. 
On returning to 190°c the second part of the curve was obtained. It will 
be seen that, in spite of the experimental errors, a time of about 10 min. 
appears to have elapsed before the initial ordering rate recovered. 
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The production of excess vacancies by cold working (Seitz 1952) can be 
expected to produce ordering at low temperatures. Figure 3.5 (a) shows 
the effect, at 100°C, of a series of light deformations performed on a partially 
ordered specimen initially quenched from 380°c. The first 500 min. 
allowed a decay of quenched in vacancies. At the end of this period a 
light deformation, at room temperature (for details of which see bd) 
produced an increase in p followed by further ordering and vacancy decay. 
Similar deformations at 500 min. intervals produced similar effects, 
although the ordering effect gradually diminished. The same experiment 
performed on a specimen initially quenched from 450°c produced no 
indication of ordering (curve (b)). This is consistent with the absence of 
domains of long range order noted above. 
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The isothermal decay at 160°c and 190°c The growth to equilibrium at 190°c 
of excess vacancies introduced by a of an initially too small vacancy 
bombardment of 1-5 x 1016 y cm~. concentration. 


Much of the data for the experiments described in § 4 e¢ seq. was derived 
from the decay of excess vacancies at 130°c. This temperature was a 
convenient one practically and had the additional advantage that ordering 
due to the thermal equilibrium vacancy concentration was barely detect- 
able. It was found that decay times were typically of the order of 10 hours 
and that the total drop in resistivity 4p, produced by an initial excess 
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vacancy concentration, was a reproducible quantity. As might be 
expected, the value of 4p, resulting from a constant bombardment was 
found to depend on the resistivity and, therefore, the degree of order. 
Figure 3.6 shows a plot of 4p, against p indicating that the ordering effect 
decreases as the equilibrium degree of order is approached. The igen 
ordered specimen used was heated between bombardments at 230°C to 
reduce p. It was noticed that, within experimental errors, the shape of 
the isotherms and hence the vacancy decay time was independent of p. 
Examples of these isotherms appear in later sections. 
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The effect on ordering at 100°c of a series of light deformations of (a) a partially 
ordered specimen, and (b) a disordered specimen, 


In order to carry out a more detailed investigation of the properties of 
the vacancies the following assumptions were made: The vacancies are 
visualized as moving through the crystal sometimes putting wrong atoms 
right and sometimes putting right atoms wrong. When the degree of 
order S is below the equilibrium value S, more are put right than wrong, 
and conversely if S>S,. The ordering rate dW/dt (where W is the 
concentration of wrong atoms) is expected to be proportional to the 
vacancy concentration v and the hopping frequency 1, and also to be 
dependent on S. The rate of change of resistivity at constant tem- 
perature is therefore taken to be 


dp dpdW_ dp ath) ; 
i aW di = Ty Unda(S)= ry Lo( W), meee tec 
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in which dp/dW and f,(S) have been combined into F,(W). On the 
usual simple models v,~v exp (—ZH,,/k7') and the equilibrium value of 
v~exp (—EH,/kT) where v is the atomic vibration frequency and E,, E, 
are, respectively, the activation energies for migration and formation of a 
vacancy. 


Fig. 3.6 
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The ordering effect at 130°c due to bombardment with 1-8 x 10'* y cm~? as a 
function of resistivity. 


It was shown in a previous paper (P1) that the ordering at low tempera- 
tures produced by an excess vacancy concentration is characterized by an 
activation energy of ~0-9ev. This is now identified as the approximate 
value of H,, (a further discussion on this point occurs in § 5). 


§ 4. H 
From eqn. (3.1), the initial ordering rate at 130°c after quenching will 
be proportional to the frozen in vacancy concentration v and the function, 
F(W). If the latter can be eliminated it should be possible, by quenching, 
to obtain the temperature dependence of v and hence #;. This was done 


544 R. A. Dugdale on some : 


in the following manner: A partially ordered specimen was queceleg 
from a series of temperatures. After each quench an isotherm at 130°C 
was taken, a fixed bombardment with y-rays given and a further isotherm 
obtained. From the quench isotherms a set of values of the initial 
gradients g, were obtained, and from the bombardment isotherms a set 7,. 
Plotting the latter as a function of resistivity allowed values of 7, at the 
beginning of the quench isotherms to be deduced by small extrapolations. 
Making the assumption that the bombardment induced vacancy con- 
centration was independent of the degree of order then gave 


Ya _“t _ const. exp (—E,/kT) Re eee 
Jn Vy 
from which H i, WaS determined. 
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Isotherms at 130°: (a) after quenching from 270-4°c, (b) after a subsequent 
bombardment with 1-4 x 1016 y em, 


In the actual experiment a partially ordered specimen was maintained 
at four temperatures in the range 250°c to 300°c for 17 hours and then 
water quenched. By way of example, the isotherms obtained after 
quenching from 270-4°c are shown in fig. 4.1 (the bombardment given was 
1-410" y-cm~*). It will be noticed that the shape of the curve is 
different in the two cases, the vacancies decaying more slowly after 
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quenching. This necessitated long heating at 130°c (not shown) before 
bombarding. On the other hand, it was found that all the isotherms of 
one type could be superimposed within the experimental error by nor- 
malizing the ordinates. This observation (discussed more fully in § 6) 
was of use since it allowed the values of g to be determined with a con- 
sistent rather than a variable error. The results obtained are given in 
the table. The final plot of log jo(g,/7,) against 10/7’ (using the corrected g,) 
is shown in fig. 4.2; a good straight line results giving H,=1-04 ev. 


: ba LO yg po LO" gy 
(°K) (4 2 cm) (4 2 cm min?) (je 2 cm) (4 2 em min-t) 
573-2 4-831 4-4 4-777 2-6 
557-6 4-611 2-0 4-587 2-0 
543-6 4-429 0-85 4-419 1-7 
529-7 4-310 0-42 4-306 1-4 
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The plot also gives 9,/g,=Vq/¥,=1 at T=—560°x. This aoe check 
on i, to be made, for, using the value just quoted, v,~5 x 107! at this 
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temperature, while, from the known y flux and data given inthe appendix, 
v,, is estimated to be ~2x 10~ (this may actually be an over estimate as 
shown in § 7.) 

A further rough check on Z, may be obtained in the following way : 
According to eqn. (3.1) the thermal ordering rate at high temperatures 
will be ae 

dt a exp(— (Ey+ Em) /kL) Fy( W). 


Now it seems reasonable to suppose that the temperature dependence of 
the ordering rate lies mainly in F(W) near to equilibrium and mainly in 
the exponential term when a long way from equilibrium. The ratio of the 
ordering rates at two temperatures 7', and 7’, for the same value of W will be 


dp dp FW) | fed 
(3), #).-<Fmee[- ( E+E) (rr ) | 


from which a quantity Q may be obtained where 


Ss a ON ce 7 dp\ |(dp om 


Putting in the ie a on we ee 6 should tend, therefore, towards 
~1-9 ev when a long way from equilibrium if the two temperatures are 
not very different. 
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This has been tested in the following manner: A partially ordered 
specimen was allowed to relax towards a higher degree of order at 190°c 
and 210°c alternately, the gradients being measured at the temperature 
change points. By making the isotherms at each temperature sufficiently 
long the small effects due to the time required for the new equilibrium 
vacancy concentration to be establiished could be neglected. Inter- 

_ mediate heating periods were given at 230°c in order that a wide range in 
p could be covered in a reasonable time. Values of Q, so obtained, are 
plotted as a function of p in fig. 4.3. The expected trend is realized, the 
highest value of Q being about 1-7ev. At the lower end, in the vicinity 


of p~4.5 wQ cm, the results show a considerable scatter. This is discussed 
in § 5. 


§5. THe Vacancy ORDERING MECHANISM 


The scatter in Q near equilibrium, just mentioned, initiated a more 
detailed study of the approach to equilibrium, this time at 230°c and 
250°C, at which temperatures the ordering rates are conveniently higher. 
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The approach to equilibrium at 250°C with alternate heatings at 230°c. The 
dashed markers on curves (a) and (b) indicate temperature change points. 
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The resistivity of a partially ordered specimen was followed as a function 
of time at these temperatures alternately. The behaviour in the early 
stages is shown inset in fig. 5.1. A measurement of Q at p~5-2 pe cm. 
gives a value of 1-4ev in rough agreement with the results described 
above. 

The interesting behaviour occurred as equilibrium at 250°C was ap- 
proached, as illustrated by curve (a) of the main graph (which starts with 
the tail of the inset curve). It will be seen that, at the 230°C to 250°C 
change points, not only does the 250°c curve have a smaller gradient than 
the 230°c curve but, closer to equilibrium, this actually becomes positive 
before again going negative. 


(a) Cb) 


The perfectly ordered structure of Cu, Au: (a) the unit cell, (b) the (111) plane. 


To explain the phenomenon it is suggested that at least two types of 
wrong atom configurations were present, one of which was able to relax 
towards equilibrium more quickly than the other. These might be called, 
respectively, easy and hard types of disorder. At the 35 hour point, for 
example, the concentration of the easy type has been brought below the 
250°C equilibrium value by heating at 230°c although this was not the 
case for the hard type. On returning to 250°c the easy type then relaxed 
upwards masking, for a while, the downward relaxation of the hard type. 

The behaviour when very close to equilibrium at 250°c is shown by 
curve (b). This was obtained immediately after (a) but is returned to the 
origin for the sake of compactness. At 230°C relaxation of both types 
occurred with the result that on returning to 250°c both relaxed upwards 
and no peak occurred. 
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To test this explanation a heat treatment was devised which, it was 
hoped, would put the easy type above the 250°c equilibrium value and 
the hard type below, i.e. the reverse of the situation which had existed so 
far. On subsequent heating at 250°c the resistivity would then be expec- 
ted to pass through a minimum. The heat treatment given consisted of 
90 hours at 230°C, producing the drop shown dashed at (c), and 20 min. at 
270°c, producing the increase shown dashed at (d). The aim was, first, to 
bring both below equilibrium at 250°c, and then to raise only the easy 
type above. On subsequent heating at 250°c the expected behaviour 
occurred as shown by curve (e). The effect was found to be quite sensitive 
to the time at 270°c, 10 min. being too short and 30 min. too long. This 
test also rules out the possibility of domain growth as the origin of the 
slower relaxation since the equilibrium domain size which would be 
implied has never been found. 

Thus the scatter on the values of Q at p~4:5 2 cm in the previous 
experiment are accounted for. The specimen must have been sufficiently 
near to equilibrium at 210°c for these effects to be significant. 

A possible explanation of the easy and hard disorder is as follows : 
Consider the alloy in the perfectly ordered state as shown by the two 
diagrams of fig. 5.2. The simplest type of disorder is produced if the 
gold atom G is interchanged with the nearest neighbour copper atom 4. 
The next simplest would be to interchange G with a next nearest neighbour 
B, and so on. Now consider how a moving vacancy might put these 
configurations right. The simplest type, the AG wrong pair, can be 
righted if a vacancy first arrives at the point G now occupied by the wrong 
copper atom. If it then interchanges with the gold atom, now at point A, 
and subsequently diffuses away (which it can do on copper sites only) 
the right configuration will result. For the wrong pair of type BG the 
process is more complicated. A vacancy has, first, to arrive at a point 
such as A and then interchange with the gold atom occupying point B. 
This now results in a wrong pair of the simpler type AG which can be 
righted either by a return of the vacancy to point G or by the arrival of a 
fresh vacancy if the first has diffused away. More complex wrong pairs 
would require correspondingly more stages in the righting process, each 
of which can be regarded as producing a wrong pair of simpler type. At 
equilibrium the concentration of the various types will be maintained by 
an equal and opposite transition between them. However, the rates of 
approach to equilibrium might well be different if the activation energies 
and other factors controlling the necessary atom vacancy interchanges 
are different. 

A similar experiment to that just described was carried out at 290°C, at 
which temperature the ordering rates are much higher, to see if a third 
relaxation process could be distinguished, but none was found. The first 
two were again evident, this time lasting for periods of ~0-2 and ~2 hours 
respectively ; the resistivity then remained substantially constant for a 
further 20 hours. It may be that the majority of the disorder is due to 
only the two simplest types of wrong pair. 
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It is well known that, for an ordered alloy such as Cu;Au, the electrical 
resistivity is not a unique function of the degree of order because of its 
dependence on both the domain size and the degree of order within the 
domains. The above observations show that even if the effect of domain 
size could be neglected p would still not be uniquely related to the degree 
oforder. This has been further demonstrated in the following experiment : 

A partially ordered specimen was prepared by two different heat treat- 
ments to about the same value of p and then bombarded after each with a 
fixed dose of y-rays (1:4 101%) em~?), and the resulting 4p,, at 130°c, 
compared. The specimen was that used for the quenching experiments 
described in § 4, and the first treatment was taken as that for the 300°c 
quench and subsequent heating at 130°c ; this resulted in p=4°78 p»Q cm 
(see the table). In this state the y bombardment produced 4p,=1-18 x 10-? 
pwQem. After heating the specimen for 2 hours at 370°c to recover the 
partially ordered state, the second treatment was given. This consisted 
of 10 hours at 230°C (producing p=4:9742 em) and 64 hours at 190°C 
(producing p=4:86u42cm). After removing excess vacancies (a small 
effect after such a treatment), the bombardment gave 4p,—0-69 x 10-? 
2 cm, which is considerably smaller. (Once again the two 130°c iso- 
therms could be superimposed ; i.e. the decay time of the excess vacancies 
was not altered.) The second treatment produced therefore, a state less 
sensitive to bombardment. It would appear that a lower concentration 
of easy disorder and higher concentration of hard disorder resulted. 

This complexity of behaviour raises the problem of what exactly is 
meant by the activation energy for migration of vacancies Z,,. According 
to the ordering mechanism here postulated there are several different 
types of hop which a vacancy may make, one of which is presumably the 
rate controlling process. Although this cannot be identified it is interes- 
ting to note that the decay time, and therefore (anticipating § 6) the number 
of such rate controlling hops before capture, is insensitive to the state of 
order. 


§ 6. THE Decay or Excess VACANCIES AT 130°C 
From eqn. (3.1) the fall in resistivity at 130°c, due to excess vacancies, 
should be 


t rt 
Ap= | +dt= | ov exp (—E,, /kT)(F W)dt. 
- 0 


Assuming that #(W) exp (—E,,/k7’) can be taken as constant if 4p is not 
too large, this may be written 


Ap=v exp (—E,, /kT)F(W) 


J 


fa 
vdt. 
0 


xs ; : : 
rhe dependence of v on ¢ is not known. One imagines, however, that the 
decay is due to the presence of vacancy capturing traps which might be 
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lattice imperfections, impurities, etc., or even interstitial atoms in the 
bombardment case. It will be Eippased therefore (following Lomer and 


Cottrell 1955) that the vacancies make an average number # of hops 
before capture such that 


A p=v exp (— E.y[kT)F(W)v ot ee OL) 


where v, is the initial vacancy concentration and ¢=nv—! exp (E,,/kT) is 
the time required for % hops. It is the aim of the present section to find 
the factors controlling 7. In the experiment about to be described the 
time ¢,). required to reach 44p at 130°c is taken as a measure of f. 

It was first thought that the domain boundaries might act as vacancy 
traps and, accordingly, an experiment to investigate the effect of domain 
size was carried out. In fig. 6.1. 4p is plotted as a function of time at 
370°C, the temperature used to grow the domains. After each period at 
370°c the specimen was cooled to room temperature fairly quickly (~5 min), 
excess vacancies removed by heating at 160°c and 130°c and 4p, obtained 
(y dose=1-5 x 101% cm~?). In some cases it was then heated at 300°c for 
2 hours and water quenched to obtain 4p,. The value of the resistivity 
after each preparation is also plotted, indicating the growth of the domains. 
Tt will be seen that there is a small dependence on domain size in the bom- 
bardment case, but hardly any in the quench case. At the end of the 
experiment, the specimen was disordered by heating at 400°C when a 
partial repeat of the process gave similar results. 


Fig. 6.1 
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Apy, Apq and p (dashed curves) as a function af time at 370°c. 


The values of #,;, for the bombardment isotherms were found to increase 
slightly from about 50min. to about 70min. although, for the quench 
isotherms, they remained constant at about 100 min. Thus, taking 
E,,=0-9ev and v=101* sec!, the values of 7% were, for bombardment 
(1,), 1-3x10° to 1:9x10° and, for quenching (%,), 2 eel U8 pment De 
observation that #,>7i, appears to be characteristic ; it was noticed in 
various exploratory experiments. An example afforded by the present 
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work comes from the quenching experiments described in § 4 in which 7, 
was 1-6 10° and %, 3:8 10°. 

The effect of domain size is therefore not very great. It was discovered 
on the other hand, that the preliminary heat treatment had a strong 
influence on 7. As described in § 2, this consisted of 2 hours at 900°c in 
argon followed by a heat treatment at 600°c. The latter was found to be 
very important in that 7 could be increased by prolonging the treatment. 

To A the effect further, a specimen was first prepared by heating 
at 900°c for 8 hours followed by cooling to room temperature in a few 
minutes (no change in weight>0-1°% occurred). The 600°C treatment 
was omitted but the ordering treatment, consisting of 5 days at 370°c 
(giving p=6-25 12cm), was given. After the usual treatment to remove 
excess vacancies, it was bombarded with 1-8 x 101° y em~* and then heated 
at 130°c to obtain a bombardmentisotherm. However, no fallin resistivity 
occurred. It was then heated at 320°c for 17 hours and quenched (giving 
p=5:35 wQ2 cm) but again, on heating at 130°c, there was no further change 
in resistivity. A quench from ~370°C (giving p=6-70 1.2 cm), which, in 
specimens given the normal preparation, would have resulted in a large 4p 

at 130°c, produced the same null effect. The values of p just quotéd 
indicate ne operation of the ordering mechanism at high temperature 
and yet no effect, due to excess vacancies, could be induced at 130°c. 
Ordering at high temperatures was verified by repeating part of the 
approach to equilibrium experiments at 250°c and 230°c described in § 5 ; 
although the ordering rates at a given p were somewhat smaller, the easy 
and hard relaxation processes were again evident. 

This behaviour is explained if the heat treatment resulted in a value of 
nso small that the lifetimes of the vacancies, and consequently the ordering 
they produced at 130°c, were negligible. This would imply that a high 
concentration of vacancy traps was induced. Their presence would still 
allow ordering at high temperatures where the equilibrium concentration 
of vacancies in solution would be maintained. __ 

The effect of the time at 600°C was next studied. The specimen, after 
being heated at 600°c for a certain time, was given a standard ordering 
treatment of 5 days at 370°c. Bombardment and quench isotherms 
(after 17 hours at 320°c) at 130°c were then taken. These, together with 
the null results just mentioned are plotted on a logarithmic time scale in 
fig. 6.2. It is apparent that 4p was influenced strongly by the 600°c 
treatment, the effect being more pronounced in the case of the quench 
isotherms. <A repeat of the 900°c treatment at the end of the ore 
again put the specimen in the ‘ null’ state. 

Values of 7%, and %, determined from these curves are plotted as a 
function of time at 600°C (log scale) in fig. 6.3. Once again %,>7, 
although both must be very small at loo0=9- The limited data suggest 
that n, saturates while 7%, is still increasing. 

It is suggested that the vacancy traps responsible for this behaviour 
are impurity atoms (not identified) which are relatively soluble at 900°c 
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and peaively insoluble at 600°c. Cooling the specimen fairly quickly 
from 900°c would then leave a supersaturated solution, the concentration 
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The effect of the duration of the 600°c treatment on the 130°c isotherms. The 
number against each curve refers to the time in hours at 600°. 
The initial resistivities of the curves, in order of increasing ¢600, were 
6:25, 6-22, 5-96, 5-90 » 2 em, for the bombardment set, and 5-35, 5-21, 
5-19, 5-20 42cm, for the quench set. 


of which could be reduced in a reasonable time (presumably by diffusion) 
by heating in the neighbourhood of 600°c. A rough estimate of the 


202 
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concentration at different stages of the experiment can be made in the 
following way : 

For a random distribution of vacancies and impurity Some the relation 
i~1/ZC would be expected, where C is the impurity concentration and 
Z is the number of neighbouring sites in which a vacancy can ‘be regarded 
as trapped. Taking Z~10 gives C~1/10% (the number of sites neigh- 
bouring a substitutional impurity is mata te 12). Applying this to the 
experimental results, the concentration was reduced to ~3x 107 by 
heating at 600°c for 342 hours. The concentration which would prevent 


Fig. 6.3 
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nm, and %, as a function of time at 600°. 

the detection of ordering after quenching from 320°c may also be estimated. 
From eqn. (6.1) 4pat. This was roughly satisfied by the quench isotherms 
for which dp~10*t), (4p in »2 cm and ¢,, in minutes). Taking the 
limit of detectability for 4p to be 5x 10-4 4Qem, this gives t,.~} min. 
and hence %,~103 and C~10-4. Such a concentration is consistent 
with the chemical analysis of the alloy. 

To account for the smaller value and possible saturation of 7, it is 
necessary to have some extra trapping sites present after bombardment 
which are not present after quenching. These might well be interstitials, 
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According to the estimate made in the appendix, the bombardment would 
put in an interstitial concentration of ~10-°. The maximum value of 7 ‘ 
reached was ~10°, requiring z in this case, to be ~10¢. This rather 
large value is interpreted to indicate that interstitials are extremely 
effective in capturing vacancies. 


-§7. Toe Besaviour or INTERSTITIALS 


Brinkman et al. (1954) concluded, from an experiment which they 
performed, that interstitials are more mobile than vacancies in Cu,Au. 
Two specimens, one partially ordered and the other disordered, were 
given identical bombardments with 9 Mev protons at —180°c which 
produced a large increase in resistivity in the first and a smaller increase in 
the second. The temperature was then raised steadily to 130°C. In 
the region —60°c to 0°c about equal annealing occurred in both, removing 
almost all the extra resistance of the initially disordered specimen. 
Approaching 130°c, a further drop occurred in the initially partially 
ordered specimen, but none in the disordered one. The low temperature 
drop was assigned to the recombination of vacancies and interstitials 
through interstitial migration, and the high temperature drop to ordering 
through vacancy migration. 

The authors pointed out that after the low temperature process was 
complete a small concentration of vacancies would remain if natural sinks 
for the interstitials existed ; due to their presence, some of the inter- 
stitials would disappear without annihilating vacancies. This obser- 
vation is here used to test their conclusions by a different type of experi- 
ment. 

Suppose there is a concentration of sinks s for interstitials and let a 
bombardment be given during the course of which 7, interstitials and v, 
vacancies are formed (i,=v,). Three interesting bombardment tem- 
peratures will be distinguished: (1) a low temperature at which both 
interstitials and vacancies are immobile, (2) a medium temperature at 
which only interstitials are mobile, (3) a high temperature at which both 
are mobile. If, after bombardment at each of the two lower temperatures, 
the specimen is raised to (3) ordering will occur. Due to the migration of 
interstitials and consequent annihilation of vacancies, the number of 
vacancies contributing to ordering will be different in each case. This 
number will now be determined as a function of dose for the three cases, 
assuming that the probabilities of capture of interstitials by sinks and 
vacancies are proportional to their respective concentrations. 

(1) At the end of the bombardment the full concentration of inter- 
stitials (1,) and vacancies (vy) will be present. During the heating to the 
high temperature the relation between the change of the one and the 
change of the other will be 

Cx) 


dv= av-+Bs di 
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where « and f are constants of proportionality. The final vacancy con- 
centration v,, when all the interstitials have gone, is obtained by inte- 
grating this expression, i.e. 
0 Vv Ss 
[a= |? Bare dy 
ty Vo av 
which gives 
20g = 10g % 
Bs VF 
(2) In this case the interstitials are mobile immediately they are formed. 


Thus the relation between the growth of the vacancy concentration and 
the formation of interstitials is 


an ae 
v=di|1—— awe 
Integrating to obtain v, 


IE iS (ae Ay 


Vp 2KV, \i? 
aE = (14 Fs ) 1s. “lal tl Ae 

(3) Here, both are mobile, the vacancies producing an amount of 
ordering dependent only on how far they get before capture. If the rate 
of production of displacements is small the vacancy concentration will 
also remain small, due to the limitation of their lifetimes by vacancy 
capturing traps (see §6). Consequently the probability of vacancy 
annihilation by moving interstitials will be small. It will be assumed 
therefore that under these conditions effectively all the vacancies con- 
tribute to ordering, i.e. 


(7.1) 


which gives 


Op==Vq bo -. Us td ye eT 


Plots of «v,/Bs against «v,/Bs, for the three cases, are shown in fig. 7.1. 
Taking ordering at 130°c to be proportional to v;, a plot of 4p, against 
dose should resemble these curves if two suitable lower temperatures are 
chosen. Such an experiment was carried out ; a partially ordered speci- 
men was bombarded at 130°c, 0°c and —196°c. The bombardment rate 
was 1-410" ycmhour!, which was small enough to satisfy the 
condition mentioned under (3) (the half time of the vacancy decay was 
~1 hour) and large enough for thermal ordering at 130°c to be negligible. 
The results obtained are plotted in fig. 7.2. The curvature of the 130°c 
plot is attributed mainly to the change in F(W) occurring during the 
course of the experiment. Figure 7.3 shows a corrected version to take 
account of this. The correction, obtained from the deviation from 
linearity at 130°C, has been applied to all three curves. 

Although the theoretical and experimental curves cannot be closely 
fitted together, which, perhaps, is not surprising in view of the elementary 
model assumed, the qualitative resemblance supports the hypothesis. 
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An estimate for s can be made by a very rough fitting, assuming the 
cross section for displacement to be 1-5 x 10-25 em? (see appendix): 
this gives s~2 x 10~%&/B. eS. 
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av,/Bs as a function of av,/Bs for the three temperatures (1), (2) and (3). 
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the three temperatures indicated. 
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According to the conclusions reached in § 6, interstititials are still present 
at430°c. Ifthisis so, some, at least, of the sinks might be better described 
as traps which anchor an interstitial without annihilating it. As pointed 
out by Lomer and Cottrell (1955), impurities might well act in this way. 
If a range of such traps exists, each with a different binding energy, s will 
be greater the lower the temperature. This may be an important factor 
in the immobilization of interstitials at low temperatures. 


§ 8. CONCLUSION 


The various experiments show that a vacancy ordering mechanism is 
a tenable hypothesis. The activation energies of formation and migration 
of vacancies, L, and £,,, are about 1-0 ev and 0-9 ev respectively. The 
interaction between vacancies and wrong atoms is complex, leading to at 
least two relaxation periods in the ordering process which can be dis- 
tinguished near equilibrium. A possible model, in terms of wrong atom 
configurations, to account for these has been given. The decay time of 
an excess vacancy concentration appears to be limited by the presence 
of vacancy capturing impurity atoms (not identified), the concentration 
of which can be varied by suitable heat treatment. Interstitials are more 
mobile than vacancies at a given temperature, but these too are probably 
trapped and can in turn contribute to the limitation of the vacancy decay 
time. 
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APPENDIX 
The Displacement of Atoms by Co y-rays 


The main interactions by which Co y-rays (mean energy~1-2 Mev) 
can displace atoms in Cu,Au appear to involve the Compton and Photo- 
electric processes. These processes lead to energetic electrons which can 
produce displacements by elastic collisions with the nucleii. The dis- 
placement cross section oc, consists, therefore, of the product of two terms 
which are : (1) the cross section for producing an electron, (2) the proba- 
bility that this will produce a displacement. An order of magnitude 
calculation, which assumes the threshold displacement energy H, to have 
the usual value of 25 ev, gives o,~10-% em?. 

Near the surface the density of displacements will, in general, be 
different to that in the interior ; electrons produced in the surface regions 
will have a chance of escape, although some compensation will occur due 
to incoming electrons produced in the surrounding material (usually mica). 
Thus, in effect, the cross section for producing displacements near the 
surface will depend on the environment. The effect was demonstrated 
with the aid of a thin specimen made from foil of thickness 0-005 cm. 
This was mounted *n a polythene container and irradiated in an 
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approximately isotropic flux of y-rays, first of all bare, and then with 
each surface covered by an additional Cu,Au foil of thickness 0-007 cm. 
Although the irradiations were both equal, the subsequent 4p, at 130°c 
was about 18°, greater after the second. The increase is attributed to a 
higher flux of incoming electrons in the presence of the extra foils, due to 
the appreciable production of photoelectrons not present in the polythene 
environment. It was estimated that about 30°/, more atoms were 
displaced (the dose was 1-4 x 1016 y cm~2). 

The extra foils in the above experiment can be considered to have 
allowed the full production of vacancies in the specimen. Advantage 
was taken of the arrangement to measure o, in terms of o, the displace- 
ment cross section for 1 Mev electrons. An electron bombardment was 
given which produced the same effect as the second y-bombardment above. 
From the ratio of the integrated fluxes it was found that o,~0,/200. 
Accepting a previously estimated value of c, to be 3x 10-*8 cm? (see P1) 
the value of cy is therefore ~1-5 x 10-?° cm?. 

Although #, has been taken ‘to be 25 ev it was found, in fact, that dis- 
placements occur with lower recoil energies. Thus 0-3 Mev electrons, 
which can only give up to 13-5 ev to the copper atoms, produced a small 
effect. Using the thin specimen mentioned above, the ratio 9.3/0, was 
found to be ~10-3 (but this is an underestimate due to the small range of 
_~0-:006 cm of 0-3 Mev electrons in Cu,Au). It may well be, as pointed 
out by Huntington (1954), that #, has not a sharp value, the probability 
of displacement depending on the recoil direction with respect to the 
lattice as well as the recoil energy. 
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ABSTRACT 

Introducing a varying degree of disorder in Mott’s quasi-crystalline 
model of liquid metals, the change in the cold neutron scattering cross 
section on melting is calculated. Motts’ calculations of the change in 
electric conductivity on melting have been repeated on the basis of this 
modified model using more recent data on the latent heat of fusion. It is 
shown that the change in the neutron scattering cross section in lead and 
the change in conductivity in sodium and potassium on melting can be 
explained if we assume that the entropy of disorder on melting is nearly 
0-5R, R being the gas constant. That this value is of the right order of 
magnitude seems to be corroborated by self-diffusion data. 


§ 1. INTRODUCTION 


RecEnt theoretical work of Placzek (1954) and others (Squires 1952, 
Kothari and Singwi 1955) has been fairly successful in explaining the 
temperature dependence of the inelastic scattering cross section of coldt 
“neutrons in polycrystalline solids. The present paper is an attempt to 
explain the cold neutron scattering at the melting point of a solid. This, 
obviously, involves a knowledge of the state of matter near the melting 
point for which no satisfactory theory exists so far. We have, therefore, 
to assume some model of the liquid state near the point of fusion. 

Mott (1954) (hereafter referred to as I) in his attempt to explain the 
change in electrical conductivity of metals in going from the solid to the 
liquid state, assumed a quasi-crystalline structure of the liquid. In this 
model it is assumed that atoms continue to vibrate as in a solid, though with 
a lower characteristic frequency due to loosening of structure, about a mean 
position which moves about at random. 

On the basis of the above model Mott obtained a reasonable agreement 
between theory and experiment. In the present paper, using a similar 
model, we calculate the change in the inelastic neutron scattering cross 
section in going from the solid to the liquid state. In contrast to the 
conductivity which mainly depends upon the Debye temperature and the 


+ Communicated by the Authors. 


{ The word * cold ’ is used here to denote neutrons of wavelength greater than 
that corresponding to the Bragg-cutoff. 
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density of electron states near the top of the Fermi distribution, both of 
which alter on melting, the situation for neutron scattering is much simpler 
in that it depends only on the former. Itis difficult to calculate the density 
of states in the liquid state since the precise shape of Brillouin zones is not 
known. 

The model that we have adopted for our calculations is that of Lennard- 
Jones and Devonshire (1939), based on the assumption that in the process 
of melting a munber of new lattice positions (called f-sites) equal in number 
to the original sites (called «-sites) appear. The redistribution of the 
atoms between the old «-sites and the new f-sites involves an increase in 
the entropy. In other words, the entropy of disorder, the nature of which 
we wish to investigate, arises from the variety of ways in which the singly- 
occupied cells in the liquid can be arranged on a pseudo-lattice. This 
entropy is a complicated function of the order-parameter which in turn is a 
function of the force constants between the atoms. With our present 
state of knowledge it is difficult to calculate the exact value of the entropy 
of disorder. 

Apart from its intrinsic interest, a study of the cold neutron scattering 
at the melting point of a solid may provide a somewhat reliable method of 
estimating the magnitude of the change in the entropy of disorder in going 
from the solid to the liquid state. The present paper is-an attempt to 
estimate the magnitude of this entropy from the electron as well as the 
neutron scattering data at the melting point. Our results are of course 
based on the assumption of a quasi-crystalline model of the liquid state. 
We do not for a moment suggest that the true state of affairs at melting 
is what is given by this model. We, however, do suggest that it appears 
to be a reasonable working model. 


§ 2. MATHEMATICAL FORMULAE 


In deriving the expression for the Debye temperature @,, in the liquid 
state, we proceed in a manner similar to that of I, where @; is obtained by 
equating the chemical potentials 1, and jz, of the solid and liquid phases 
respectively, at the melting point 7’). ; 

The chemical potentials pg (7'y,) and uz(7");) are respectively (see Fowler 
and Guggenheim, Statistical Thermodynamics, p. 329) 


bs (T= — x99 +3kT y log, [1—exp (—hvs/kTy)|—ET y log, Js (Tm) 


and =p Vs, 2 . . . ° A . . i (1 a) 
+pV,—akT y, 4 : 6 “ ° . cj (1 b) 


where the subscripts L and S denote the liquid and solid phases respec- 
tively, L is the latent heat of fusion, / is the Boltzmann constant and other 
terms have their usual significance. ore 
Formulae (1a) and (16) are based on the Einstein approximation 
which is valid in our case since 7’>@. The additional term «k7'y in 
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(1 b) represents the energy corresponding to the entropy of disorder in the 
liquid state arising from the number of ways in which the singly-occupied 
cells of the liquid can be arranged on a pseudo-lattice. In fact «k repre- 
sents the entropy of disorder and is given by eqn. (1) of Appendix A, the 
derivation of which is given by Lennard-Jones and Devonshire (eqn. 
(12-B), 1939). We have here assumed « to be independent of temperature. 

On equating p(Z',) and p,(7'y,) and neglecting the small differences 
between both the third and fourth terms of (1a) and (1b), we obtain the 
following formula for @,, 

L—akT y ee 1—exp (—9,/T'y) 
SkT ye ao Leen Ce Opray 

where L=y,°— y,°, represents the latent heat of fusion. 

Equation (2) differs from eqn. (1.1) of I in that the former contains an 
extra factor «kT. It is difficult to estimate the value of « in the absence 
of any proper theory of melting. We have, therefore, used in our cal- 
culations different values of «, ranging from 0 to 1 to see which value 
gives a reasonable agreement with experiment. 

If we use Placzek’s incoherent approximation the inelastic neutron 
scattering cross section is given by the formula 


‘ __ 3(S+s) ieee 2 ie en /S ieee! 
"ru Mk, \(-7+3¢ tS @, 2 


TTR Cel aa Le 
totem) tet + 
(see eqns. (40) and (41), Kothari and Singwi 1955). In eqn. (3), S and s 
denote respectively, the coherent and the incoherent scattering cross 
sections per atom, V is the atomic mass expressed in terms of the neutron 
mass my and k, is the wave number in dimensionless units corresponding 
to the neutron wave length A and is given by 


Rh? \U2 Qa 
C= ae ees es ee 
ae (=-3) OgV2 r (4) 


We assume that the scattering cross section in the liquid state is also 
given by (3) and (4) with @, replaced by @,, and where @, is given in 
terms of @, by (2). 

It may be pointed out here that though formula (3) takes care of all the 
multi-phonon processes which are important for temperatures of interest, 
its derivation is based on the assumption that atomic vibrations remain 
harmonic. The neglect of the contribution arising from anharmonic 
effects is a simplification. Since we are interested only in the difference 
of the scattering cross sections in the liquid and solid states at the melting 
point, the neglect of anharmonic contribution may not be serious. 

In order to see whether the value of « as obtained from the neutron 
scattering data is consistent with that obtained from the change in 
electrical conductivity! in going from the solid to the liquid state, we 
have repeated Mott’s calculations, using recent values of the latent heat 


(2) 
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of. fusion and also taking into account the factor «kT'y. Without re- 


peating Mott’s arguments we shall write the expression for the ratio of 
electrical conductivities as 


(78) ie (3) [K(dB/dK)] <2 : 
Ue GR) e aes oe A 


where the subscripts L and S have the same meaning as before, K denotes 
the wave number of an electron at the top of the Fermi distribution and 
EF is the kinetic energy of such an electron. In writing eqn. (5) we have 
assumed with Mott (1934) that all other factors in the usual conductivity 
expression except those retained on the right-hand-side of (5) do not 
change appreciably on melting. Our model of the liquid state enables 
us to calculate ©, in terms of @,, which is known. It, however, does not 
enable us to calculate [K(dH/dK)], in terms of the corresponding quantity 
in the solid state, thus introducing considerable uncertainty in the cal- 
culation of y,;. For good conductors and especially for monovalent 
metals, using the same argument as given by Mott, we can write (5) with 


good approximation as ns _ Os a 
Alias Tc nei ea 


The validity of (6) is doubtful for other metals. 


§ 3. RESULTS 
Using eqns. (2), (3), (4) and (6) and different values of « we have cal- 
culated (7 s/nz)r,, and (7,—c¢g)/og x 100 for various metals, and the results 
are shown in table 2. 
Values of constants used in our calculations are given in table 1. 


Table 1. Data for Various Metals 


Dene Wave number of 


neutrons in 
dimensionless 
units (eqn. (3))t 


Bs 


Melting | Latent heat 
point | of fusion in 
Ty °K cal/gm 


Atomic 
weight 


Temp. S+s 
(solid) |(barns) 
Og °K 


Element 


207-2 600-4 ° 5-86 88 0:3963 
(8:3) 

63:54 | 1356 g 0-2191 
(8-0) 

26-97. | 932-7 f 0-1933 
| (8-0) 

24-32 | 924 cf 0:2827 

(6-25) 

22-997 Be ; 0-2939 
(8:3) 

9.096 es 0:3736 
39 pee 


+ Quantities in brackets give the corresponding wave lengths in A units. 
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It will be seen that for «0-5 the calculated values of (ns/nz)ry, for 


Na, K and Al are in good agreement with experimental values, whereas 
for Pb and Cu «0, gives a good agreement. Unfortunately, experi- 
mental data on the change of cold neutron scattering cross section at the 
melting point are lacking. However, preliminary results of Squires} on 
lead indicate a sharp increase of about 20° in the scattering cross section 
at the melting point. This would seem to indicate a value of « of about 
0-5, as can be seen from table 2. On the other hand the conductivity 
data require « to be about 0. 


§ 4. Discussion 


The discrepancy between the two values of « in the case of lead is 
perhaps due to our assumption that the value of [K(dE/dK Niw is the 
same both for the solid and liquid phases. There is at present no possible 
way to decide this point quantitatively. We may, however, remark that 
a more reliable estimate of the value of « is obtained from thermal in- 
elastic neutron scattering rather than from conductivity data, since the 
latter involve an uncertain change in [K (d#/dK Near 

The situation in the case of alkali metals is somewhat better. As we 
see from table 2, for Na and K the calculated values of (79/7 L)ry ere in 
good agreement with experimental values for «0-5. Now the free 
electron approximation is known to hold for these metals; hence the 
change in conductivity on melting can be ascribed wholly to the change 
in their Debye temperatures. Unfortunately, no neutron scattering 
data are available for these metals to decide conclusively in favour of 
“S05. 

The conductivity data of Na and K and the neutron scattering data of 
Pb seem to indicate a value of «0-5 for these metals. This would 
mean that the change in the entropy of disorder on melting for these . 
metals is nearly 4R, R being the gas constant. The extra entropy is 
undoubtedly a manifestation of disorder in the liquid state and arises 
from the variety of ways in which the singly-occupied cells in the liquid 
can be arranged on a pseudo-lattice. It might be thought that this 
entropy of disorder in the liquid state could be explained by allowing the 
possibility of more than one atom to occupy the same cell. Recent cal- 
culations of Pople (1951) have shown, however, that these states give only 
a negligible contribution to the liquid partition function. The entropy 
arising from these states is called the communal entropy which does not 
become important until the critical point of the liquid. 

We shall now use the order disorder theory of melting of Lennard- 
Jones and Devonshire (1939) to see whether the value of « as obtained 
from the conductivity and neutron scattering data is of the right order 
of magnitude. The relevant equations are given in Appendix A. 


+ We are grateful to Dr. Squires for sending to us his experimental results 
on lead before publication. Dr. Egglestaff informs us that he also observes a 
similar change. 
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In their theory the entropy of disorder and the average energy W 
required to take an atom from an «-site to a B-site are given in terms of 
an ‘order parameter’ Q. Taking the entropy of disorder to be 0-58, 
corresponding to «0-5 and using eqns. (A.1) and (A.2), we find that at 
the melting point the values of W for Pb, Na and K are respectively, 6-2, 
3-8 and 3-4 kcal per g mole. The «- and f-sites in the theory of Lennard- 
Jones and Devonshire have been introduced merely to obtain a con- 
venient representation, capable of mathematical treatment. The argu- 
ments of the theory will equally apply to the general case of the trans- 
formation of a close-packed crystalline array to the more disorderly, open- 
packed arrangement in the liquid. It, therefore, seems quite reasonable 
to associate the average energy W with a part of the activation energy for 
self-diffusion in the liquid state. 

The self-diffusion data in liquid lead at the melting point (see Appendix 
B) give for the activation energy a value of nearly 14 kcal per g mole. 
Self-diffusion in Pb is presumably determined by vacancies if the work 
of Huntington and Seitz (1942) on copper is taken as a basis since as in 
the case of copper the inner shells in lead are large. These authors have 
shown that in copper the energy of formation of a vacancy is of the order 
of the energy of migration. If we assume the same thing for lead, then 
the energy of formation of vacancies is nearly 7 kcal per g mole. This 
value is not far from 6-2 kcal per g mole as obtained by us. The agree- 
ment is more than what one would expect on the basis of our crude model. 

In Appendix B we have estimated the activation energy for self- 
diffusion in liquid Na and K to be nearly 7 keal per g mole. Huntington 
(1942) has shown that the energy required to form vacancies and that to 
form interstitial atoms are much more nearly the same for alkali metals 
than for copper. As a result it is not possible to decide between inter- 
stitial and vacancy diffusion in alkali metals. Whatever may be the real 
mechanism of self-diffusion in alkali metals, it is gratifying to see that 
values calculated by us are of the same order of magnitude as those given 
by self-diffusion data. 

We, therefore, arrive at the conclusion that the conductivity data of 
Na and K and the neutron scattering data of Pb require that on melting 
the entropy of disorder for these metals is nearly 0-5R, and that this 
value is consistent with self-diffusion data. It would be highly desirable 
to have experimental data available for cold neutron scattering cross 
section in metals near the melting point. 

In trying to explain the abrupt change in cold neutron scattering cross 
section at the melting point in lead, we have neglected the possibility of 
elastic scattering contributing in the liquid state. The neutrons used in 
Squires’ experiment were filtered through graphite, which gives a cut-off 
at 6-7 A, the mean wavelength being 8-3 A. The cut-off in lead is at 5-7 A. 
Now the observed change in volume on melting in this case is less than 5°, 
which would correspond to a change of utmost 2° in linear dimensions. 
Thus on our model, graphite filtered neutrons with a cut-off at 6-7 A would 
not give rise to Bragg reflections in liquid lead, 
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Finally we would like to add that the scattering of cold neutrons in 
solids may provide a powerful method of investigating the absolute 
number of defects produced either by temperature or by irradiation. 
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APPENDIX A 


In the theory of Lennard-Jones and Devonshire (1939) the entropy 
of disorder is given by 


S=—2 (1—@Q) log (1—Q@)—2Q log Q, . . . . (A.1) 
where the order parameter Q is defined by Q@=N,/N, NV, being the number 
of atoms on «-sites and N is the total number of atoms. And the average 
energy W required to take an atom from an «-site to a f-site is 


W=4kT tanh-1(2Q—1). . . . . . (A.2) 


APPENDIX 3B 
The coefficient of self-diffusion is given by the usual relation 
D=D, exp (—E/RT) Rte tk tsk) 


where £ is the activation energy per g mole. In the case of (solid) lead # 
is 28kcal and D, is 5-1 cm?/sec. Using the measured value of 
D=2-55 x 10~-* cm?/sec in liquid lead (Groh and Von Hevesy 1920, Ann. 
der Physik, 63, 85) near the melting point, and assuming that the value of 
_D, does not change appreciably in going from the solid to the liquid state, 
we get the energy of activation # in the liquid state to be 14-5 kcal. 

In the case of alkali metals Na and K the activation energy H# in the 
solid state from self-diffusion data is nearly 10 kcal and the value of D, 
is 0-32 cm2/sec (Norberg and Slichter 1951, Phys. kev., 83, 1074, who quote 
private communication from Nachtrieb). Since in most metals in the 
liquid state the coefficient of diffusion D is of the order of 10~° cm*/sec, 
we assume this value for Na and K and as before take the value of Dy 
corresponding to that in the solid state to calculate the value of # in the 
liquid state. The calculation yields 


Element Na K 
E 7:6 6-9 
(keal/g mole) 


S Va vi INE =6 7p 
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LVII. Creep in Metal Crystals at Very Low Temperatures 


By N. F. Mott 
Cavendish Laboratory, Cambridget 


[Received December 20, 1955] 


SUMMARY 


A theory is given of creep in metallic single crystals at very low 
temperatures, which is based on the hypothesis that dislocations can 
pass through barriers, such as those provided by dislocations crossing the 
slip plane, owing to the quantum-mechanical tunnel effect. A detailed 
theory would depend on a knowledge of the atomic arrangements at the 
centre of a dislocation, which in this case does not exist, but a rough 
calculation shows that creep should be independent of temperature 
below about 10°K in cadmium. 


In a recent paper in this journal Glen (1956) has described measurements 
of creep in single crystals of cadmium at the temperature of liquid helium, 
and has confirmed former observations of Meissner, Polanyi and Schmid 
(1930) that the effect exists. The purpose of this note is to attempt 
a theory of creep at low temperatures, based on the idea that a dislocation 
held up by a potential barrier, such as that which is formed when a screw 
dislocation cuts the plane in which it moves, can pass through it owing 
to the quantum-mechanical tunnel effect. 

The model is that used to explain logarithmic creep (see, for example, 
Mott 1953). An (edge) dislocation is held up by a series of obstacles 
S,. S,, 83; which may for instance be screw dislocations cutting the plane 
of the paper, so that the dislocation cannot move forwards without 
forming a jog. It is assumed that the stress o acting on the dislocation 
is nearly as great as the stress oy required to move it without the help 
of temperature. It is usually supposed (cf. Mott 1953), Cottrell 1952) 
that the activation energy to form a jog, and thus to pull a dislocation — 
through a barrier, is then 

J(l—ofog)s~ + oh nace, 1g 


where J is the jog energy. If one supposes, however, that the energy 
of the jog is a smooth function of displacement x along a line perpendicular 
to the dislocation, as in fig. 2, a better approximation would be 


Wool og aon ovate eee 


and the width of the barrier 


w=a(1—o/o,)'. ve pte ok) sk Aen red 
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(a) Showing two stages in the movement of a dislocation line. 
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(a) is the potential energy of a jog as a function of displacement of the 
dislocation. 


(b) is the same with the term oblx due to the stress added to it. 
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Here W, is of order a few electron volts and a of the order of the inter- 
atomic distance. 

Expressions of this general form were first given by Mott and Nabarro 
(1948) for a different model, namely for dislocations held up by the stress 
field round a precipitate, but it seems to have escaped notice that they 
are of more general applicability. We doubt whether the change will 
greatly effect predictions already made except when o is very near to 
o), and hence at very low temperatures. We should under these 
conditions expect ‘logarithmic’ creep to satisfy the law (Mott 1953, 
p. 752) 
(strain)?/?=const. k7' In (yt), 


and hence the flow stress to depend on temperature according to a law 
of the type 
Mey o=0y—AT32. 


In the problem to be considered here the possibility of the tunnel 
effect arises through the fact that, when o approaches og, the barrier 
becomes low and narrow. A moving dislocation is thought to have 
effective mass of order M, the atomic mass, per atomic plane cut by it. 
Tunnel effect will only be possible in processes in which a small length 
of dislocation moves. We thus suppose that, in creep by tunnel effect, 
a jog is formed by a process such as that shown in fig. 1 (b); the dislocation 
moves from the form AS,B to AS’,B, S,B being a length equal to na so 
that the moving dislocation cuts 2n planes. The force acting at S’, 
is smaller than that acting at S,. We suppose that the effective mean 
stress between S, and S’,, points distant w apart, is 


Oop (l—w/na)o. «! ae eee 
We thus have an increased barrier height compared with (1 6) 
W= W)(1—o.g/og)* Ree 
and width compared with (2) 
w=a(1—ogg/o9) 1/2. eh 4 eee at 


The latter equation with (3) gives a quadratic for og; we have from 


(3) and (5) ‘i , Ba are 
ely \ ae es Be Oo. 
( e 2) een (: ca ; 


To a rough approximation suppose » to be large enough to permit an 
expansion Of og in ascending powers of 1/n; we then find for the factor 
that occurs in (4) and (5) 
Oo ] 
jaa — liz 
To Rp rn 


where n=1l—da/op. 
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We now take the chance per unit time, to a sufficient approximation, 
that in fig. 1 (6) the dislocation moves by tunnel effect from the con- 
figuration AS,B to AS’,B to be 


P 4n MW 1/2 
v exp —2( Re wh eC eee a ee ae 0) 


this being the chance for a ‘particle’ of mass 2nM to penetrate a 
(rectangular) barrier of height W and breadth w. : 

v is an atomic frequency of order 1012 sec-! and whose exact value 
need not concern us. The dislocation is treated as a dynamical system 
with many degrees of freedom; penetration will occur for the mode of 
motion for which the value of x makes (6) a maximum. The expression 
(6) is a maximum when v is chosen to give a minimum value of 


(nW)*w 


n2(1— eq /g)/4= n(n 5/4, 


and hence of 


Since 7 is small it is clear that this minimum value is reached when 
n=1/n/?, so that (6) becomes 


yexp| —2 iG an} 


where « is a constant of order unity. 
The creep strain « will be given by a law of the type : 


e=A In yt 
where the constant A is independent of the temperature and given by 


A=const./a,/(2M W,/h?) 


and the constant is of order unity. 

We shall compare the creep rate with that due to thermal activation. 
According to the analysis given above, the chance per unit time for 
a thermally activated movement will be 


expo eign il ts.) 4.0 Meet. ce © (7) 
Creep by tunnel effect will thus be more rapid than by thermal activation 
if 
2MW, Wie ara 
| a 4 (=F) rein 


A reasonable assumption is that 2mW, a?/h?~1, where m is the mass of 
the electron ; this would imply that W, is about lev. W, should be of 
the order of the jog energy, but without a detailed knowledge of the 
curve shown in fig. 2 we cannot evaluate it directly. The critical tem- 
perature above which thermal activation becomes important is then given 


by kT |W yh (rn Mh. 
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Since the creep rate is observable, the quantity (7) should be of the order 
of the inverse of a few seconds, so that W n*/?/k7' will be of the order 30. 
Some re-arrangement then gives, eliminating 7 


kT m \3/4 (1 \3/2 
(GY Grom 
0 a 


If m/M~0-5x 10-> and W,~l ev, this gives temperatures of the order 
of a degree Kelvin. The various approximations introduced clearly 
make this result uncertain by several multiples of two; it at any rate 
predicts creep independent of temperature in the liquid helium range. 
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LVIII. The Lifetime of the 200 kev Excited State of 72 


By C. M. P. Jounson 
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EXPERIMENTS to determine the magnetic moment of the 200 kev excited 
state of '°F (see Phillips and Jones 1956 and Lehmann, Leveque and Fiehrer 
1955, and Treacy 1955) have made it necessary to know the lifetime of the 
state more accurately than previously determined (Jones, Phillips 
Johnson and Wilkinson 1954, and Thirion, Barnes and Lauritsen 1954). 

The experimental arrangement used was similar to that of Bell, 
Graham and Petch (1952) as shown in fig. 1. The counters were sodium 
iodide crystals 2 in. x 2 in. x 1 in. in size, optically coupled to EMI 6262 
photomultipliers. The collectors of these were connected to limiters for 
the fast coincidence circuit, and proportional pulses were taken from the 
twelfth dynodes. The resolving time of the fast coincidence unit was 
45 millimicroseconds in one series of runs, and 80 millimicroseconds in 
another. The proportional pulses from the counters were fed to single 
channel kick sorters, set for the counter A to accept the photopeak of a 
1-35 Mev gamma ray, and the counter B the photopeak of a 200 kev 
gamma ray. The outputs of the single channel kicksorters, and of the 
fast coincidence unit, were connected to a triple coincidence circuit of 
resolving time about one microsecond. ‘This ensures that each counter 
counts only one gamma ray. 

The sources of 19F used were discs of teflon irradiated for one minute with 
fast neutrons from the cyclotron, giving 190 which decays to 1°F with a 
30 second half life. The discs were placed between the two crystals after 
irradiation, and the coincidence rate recorded for different lengths of 
delay cable inserted between counter A and the coincidence unit. Coin- 
cidence counts were taken for the same number of single counts from the 
low energy counter, so that the number of decay processes in each run was 
the same. The runs were started when the activity of the source had 
reached a standard value, so that the number of random coincidences was 
the same for each value of decay, and any dead time effects in the singles 
scaler were eliminated. The number of random coincidences was deter- 
mined by putting the delay between counter B and the coincidence unit, 
and was found to be 4 for each run, so that the ratio of real to random 
coincidences was never less than sixteen, even for the largest value of 


delay. 
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The effect of attenuation in the delay cable was investigated by plotting 
prompt resolution curves using a ?2Na source first with no additional 
delay, and then with 40 metres of cable in each counter lead. A small 
decrease in the effective resolving time was observed, hence a fairly long 
resolving time was used so that the effect of this on the lifetime will be 
very small. This was confirmed by the fact that there was no significant 
difference in lifetime from the runs with the two values of resolving time. 

The delay cable was calibrated by measuring the resonant frequencies 
of the three lengths of cable used with an accurate crystal calibrated 
wavemeter. ‘This introduced a correction of 2°% to the lifetime since the 
cables were cut to length assuming the manufacturer’s value of dielectric 
constant to be correct. A semilogarithmic plot of the sum of all the 
runs, with randoms subtracted but not corrected for cable error, is shown 
in fig. 2. The coincidence rate for zero delay was not used as there may be 
spurious prompt coincidences due to the decay of 1®N or detection of 
beta—gamma coincidences. 

The slope of the line was calculated by the method of least squares to be 
(0-812-+0-0087) x 10% sect, (probable error) giving, after correction for 
cable error.and allowing for the small uncertainty due to attenuation, a 
value of the mean life of (1-25+0-025) x 10-* sec. This agrees well with 
other determintaions, (Fiehrer ef al., private communication, Thirion 
et al. 1954) and with the matrix element calculated from the shell model 
by Elliott and Flowers (1955). 


The author is grateful to the Department of Scientific and Industrial 
Research for a maintenance grant. 


REFERENCES 


Bett, R. E., Grauam, R. L., and Petron, H. E., 1952, Canad. J. Phys., 30, 35. 

Exriorr, J. P., and Frowers, B. H., 1955, Proc. Roy. Soc. A, 229, 536. 

Jonzs, G. A., Puitiipes, W. R., Jonnson, C. M. P., and WILKINSON, IBY, Jake, 
1954, Phys. Rev., 96, 547. : 

Lenmann, P., Leveque, A., and Frenrer, M., 1955, Comptes Rendus, 241, 700. 

Pumps, W. R., and Jonss, G. A., 1956, Phil. Mag., Sune. 

TurrIon, J., Barnes, C. A., and Lauritsen, C. C., 1954, Phys. Rev., 94, 1076. 

Treacy, P., 1955, Nature, Lond., 176, 923. 


LIX. The Magnetic Moment of the 200 kev Excited State of °F 


By W. R. Puitires and G. A. JONES 
Javendish Laboratory, Cambridgey 


[Received December 22, 1955] 


SUMMARY 


We have measured the magnetic moment of the 200 kev excited state 
of 19F by observing the influence of an external magnetic field on the 
directional properties of a yy cascade through this state. The cascade 
follows the B-decay of 180 produced by neutron bombardment of HF. 
The result, p=+3-50+0-24n.m., assumes that no time-dependent 
quadrupole interaction is present in the source ; the effect of the possible 
presence of such an interaction is discussed. 


~§ 1. InrRopUCTION i 


THe interaction of the magnetic moment of an excited nucleus with an 
external magnetic field may alter the angular correlation of radiation 
proceeding through the excited state. Observation of the change may 
then be used to determine the magnetic moment of the excited nucleus 
and such measurements have been performed, e.g. Aeppli et al. (1952), 
Raboy and Krohn (1954). For the change in correlation to be observable 
with present techniques the lifetime of the excited state must lie within 
‘the region 10-® to 10-* sec. 

We have measured the magnetic moment of the 200 kev excited state of 
19F by observing the change in correlation with magnetic field of the y—y 
cascade proceeding through this level after the B-decay of 1°O (Jones 
et al. (1954). The mean lifetime of this state has been determined as 
(1-25 -+0-03) x 10-7 see (Johnson 1956), and its spin is 5/2. 

Liquid sources of hydrofluoric acid were irradiated with neutrons from 
the Cavendish cyclotron to produce the 1°O. The lifetime of 30 sec 
allowed removal from the vicinity of the cyclotron to make the observa- 
tions. The measurements have been analysed on the assumption that 
there is no perturbing influence in the intermediate state other than that 
of the external magnetic field, but the effect of such interactions it 
discussed. 

During the course of this work, similar measurements, utilizing different 
reactions involving this state, have been carried out at Saclay (Lehmann 


et al. 1955) and at Canberra (Treacy 1955), with results in excellent agree- 
ment with our own. 
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§2. THEORY 


The interactions of magnetic and electric fields on the magnetic and 
electric moments of the nucleus induce transitions between the magnetic 
substates in the intermediate state, defined with respect to an arbitrary 
axis. The resulting time-dependence of the population of the substates 
implies a time-dependent angular correlation of the emitted radiation. 
In general if W,(@)=1+ 2, A, P,,(0) is the unperturbed correlation, then the 
perturbed correlation is given by W(#)=, A ,G, P,(0) where the coefficients 
G,. of magnitude less than unity, are given for several cases of interest by 
Abragam and Pound (1953). For the second y-ray emitted at time ¢ after 
the first, the delayed correlation is 


Weve te SA GP (6).-5. .-. . « 2-101) 


In the case that the only interaction present is that of the magnetic 
moment j with an external magnetic field H it may be shown that 


Wee ja Vireo’ Wabtewt) . . . . . « (2) 


if coincidences are observed ‘one way only’ (i.e. one detector has zero 
efficiency for one of the y-rays) as in the present experiment. 7 is the mean 
lifetime of the intermediate state and w is the Larmor frequency of pre- 
cession, 1.e., o=wH/Ih=gu,H/h where g is the gyromagnetic ratio and J 
the spin of the excited state. 

tf the coincidence resolving time is very much greater than the lifetime 
of the state than the average correlation is observed : 


W(6, H)= | W(t) dt. 


0 
For W,(@)=1-+A, cos? 6 this leads to 


Ww (180°, Hf) ie 1462+ 4w?7? m 
A 
where bo= gt. 


Thus a series of measurements of this ratio at different values of H gives 
||, whilst observation of W(0, H) for a particular H gives the sign of 9. 
Because of poor yield, this was the method used to determine g, rather | 
than the more elegant observation of W , (e.g. Steffen and Zobel 1955). 


§3.-Decay or 120 


The features of the decay scheme of 1°O pertinent to the discussion 
are shown in fig. 1. The 200 kev state is known to be (5/2+-) (Jones et al. 
1954, Sherr ef al. 1954). We assign the characteristics (3/2+4-) to the 
1:59 mev state from measurements, with liquid HF source, of the unper- 
turbed y,-y, correlation, which is adequately represented by the function 
W,(0)=1—(0-272-+0-027) cos? 6, after correction for geometry. Of the 
three possibilities, (3/2+), (5/2++) and (7/2--), consistent with the decay 
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scheme, (3/2) gives exactly this correlation for a pure MI transition 
for y,; (7/2+-) can be made to fit the data only by assuming a 3°, 
admixture of E2 in y,, i.e. about’ 300 times more than for single particle 
transitions. The most recent survey of radiative transitions in light nuclei 
(Wilkinson 1956) indicates that such an admixture, though not very 
probable, cannot be ruled out. We prefer the assignment (3/2+-) because 
shell model calculations on this nucleus (Elliott and Flowers 1955) 
predict a (3/2++) state near 1-6 Mev with a strong M1 transition to the 
200 kev state. The spin assignment does not affect the interpretation 
of the experiment, but is of interest in assessing possible perturbing 
effects in the source. 


§ 4. CHOICE OF SUITABLE SOURCE 


In the previous note on the !°O decay it was observed that, using teflon 
sources the observed anisotropy was two-thirds the value expected for 
y, a pure M1 transition to a (3/2+-) state. Further measurements showed 
that in fact the anisotropy was initially the same as that observed with 
HF sources but changed on neutron bombardment of the teflon, decreasing 
to about one-fifth the initial value. The structural disorder caused by 
the neutrons thus introduces a large perturbing effect, and the final 
anisotropy is consistent with the © hard core’ correlation remaining in 
the source when a static electric quadrupole interation is present. 

Using 40%, by weight hydrofluoric acid as source the only expected 
interaction is that of the electric quadrupole moment with randomly 
fluctuating electric field gradients. For such an interaction Abragam 
and Pound (1953) show that the coefficients G,(f) m eqn. (1) have the 
values exp (—A,#). A, is proportional to the correlation time 7, of the 
HF solution, which is a measure of the time the environments of a 
molecule in the liquid are in a given configuration. 7, is roughly propor- 
tional to the viscosity of the liquid and as 7 ~ 1-5 cpoise in our case, we 
expect A, to be very small. 


§5. EXPERIMENTAL MEetTHop 

The y-detectors were Nal (Tl) crystals, 3-7 em dia. x 3-7 em long for 
the 1:37 Mev radiation and 2cmx2em x1 cm for that of 200 kev. To 
assist magnetic screening, 6 in. long perspex light-pipes removed the 
accompanying photomultipliers type EMI 6262 from near the magnet. 
Perspex absorbers, to remove f-radiation, were placed over the detecting 
crystals, and, to prevent the possibility of spurious coincidences by 
y-scattering from one crystal to the other, the high energy detector was 
screened, front and sides, with 4 mm of lead. 

The electromagnet was designed with the object of minimizing the 
occurrence and effect of spurious coincidences arising from the large 
flux of B-rays stemming from the source and in Coincidence with Vi 
and y2. ‘To this end the pole pieces were long and tapered (3 in. dia. 
to 1 in. dia. over 6 in.), and the yoke and coils were removed far from the 


¢ 


+595 nev. 


1-45 nev 
ad 37 mev 
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region of detection. In addition the pole-faces were capped with graphite 
discs, 1 em thick, to reduce the intensity of bremsstrahlung ; the remain- 
ing gap of 2cm was sufficient to accommodate the large sources, in 
Papece polythene containers of inner dimensions 1 one dia. x 1-5 em. 

e variation of H over the ga ay vn SS 
than 1°, over the dimensions of ie pan Ee eae 
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Decay scheme of 1°O, relevant to the Schematic diagram of electronic apparatus 
experiment. used. 


A block diagram of the electronic apparatus is shown in fig. 2. The 
detecting ‘ windows’ of the single channel kicksorters were set respec- 
tively to straddle the photopeak of y, and to cover the range 600-1400 kev 
of the spectrum of y,. The efficiencies of y-detection were maintained 
constant by frequent monitoring of the relative positions of * window ’ 
and peak using a multi-channel kicksorter in conjunction with a gate 
circuit. The mean value of the coincidence resolving time, determined 
throughout the experiment, was 27),=0-7-+-0-01 psec. 

It was found that the reaction 19F(n, «)!®N produced an intense source 
of 7-sec 16N, but that a waiting period of 1 minute between irradiation 
and measurement eliminated the possiblity of significant effects from this 
decay since y-cascades are not produced (see Ajzenberg and Lauritsen 
1955) for y-transitions in 160. 

The total number of coincidences observed divided by the total number 
of, say y, single counts at the various angles is a measure of the correlation 
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function provided that the efficiency of the y,-counter does not change. 
In our case the y,-counter was moveable and the efficiency of detection 
of y, was kept constant by monitoring with the multichannel kicksorter. 
If N, is the number of coincidences observed, after subtraction of randoms, 
for a total of N counts in the low energy detector, after subtracting the 
background, then R(#) is a measure of the correlation at field H where 


W(180°, H N N, 
R= = (st) | (a) 
W (90°, #) N iso | \N J oor 
The random coincidences were estimated from the single counting rates 
and the resolving time, allowing for the decay of the source ; they never 


exceeded 5°, of the total number of coincidences. The values of R(H) 
for different values of H were measured in addition to the zero-field ratio. 


$6. RESULTS 


Before the results obtained can be analysed in terms of eqn. (2) the 
following corrections must be considered :— 


(1) Geometry corrections ; eqn. (2) refers to a point to point correlation. 
Krohn and Raboy (1955) have shown that the geometry correction is 
independent of field, but only for circular detectors of constant efficiency 
over their surface. It is not expected that the small departure from these 
conditions in our case will seriously affect this conclusion, so the geometry 
correction is simple. 


(2) The effect of finite resolving time ; eqn. (2) is derived for infinite 
resolving time. We have estimated that, for values of H of interest, 
the correction required is never more than }°,, so we have neglected it, 
but included its effects in determing the probable error of our result. 

(3) The effect of spurious coincidences produced by the presence of 
the magnet. The zero-field ratio R(0O) was greater than the value obser- 
ved on removal of the pole-pieces by an amount corresponding to a 10% 
isotropic contribution of coincidences from secondary effects in the 
pole-pieces, and this istoropic contribution must be subtracted out from 
the experimentally observed ratios R(H). Since all the ratios R(H) 
contain information on the ratio R(O), we have attempted to arrive at 
a correction based on all the observations rather than merely on the ratio 
of R(O) with and without the pole-pieces. ; 

The method of analysis was as follows. A value was assigned to R(0), in 
a range consistent with the experimental measurement. and the ‘magnet 
correction’ therefore determined. A least squares fit for the corrected 
(fH) to eqn. (3) was made in the parameters b, and wr and this was 
compared at H—0 with the corrected R(0), namely the experimental 
measurement without pole-pieces. Over a range of values of R(0), the 
least squares fit was consistent at H=0 with the corrected R(O) and the 
curve corresponding to the centre of this range was taken as the best 
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possible fit, with the appropriate w7 as the best experimental determina- 
tion. This curve is shown in fig. 3, together with the corercted R(/) ; the 
errors on the points are probable errors based on counting statistics. The 
value of gz calculated from this curve is g7=1-74 107%, which, for 
+= 1-25 x 10-*, gives | g|==1-40+0-09. 

The sign of g was determined by observing W(6, H) for H=700 gauss 
The correlation is shown in fig. 4. The direction of precession is that 
corresponding to positive ju, so, since the spin J=5/2, we have 


w=+3-50+0-24 n.m. 


The errors quoted are probable errors and include errors arising from 
uncertainty in b,. spread of experimental points about the curve of best 
fit, and uncertainty in 7 and H. 


$7. Discussion 

The analysis of the results in terms of eqn. (2) assumes that there is no 
time-dependent quadrupole interaction present in the source. Such 
interaction would decrease the anisotropy of the correlation, so an E2 
component would be needed in y, to offset this and the agreement between 
measurement and a pure radiation pattern reduced merely to chance. We 
have calculated the effect of an interaction, of strength A,= 1-83 x 10® see~}, 
such that the true anisotropy is 20°, greater than that measured at H=0, 
and find that g would be increased by 25°%, while the necessary introduc- 
tion of 4% of E2 in y, is probably not unreasonable. Thus the effect of 
quadrupole interaction in the source is obviously of importance. So far, 
because of poor yield, we have not attempted delayed correlation experi- 
ments which would indicate the magnitude of such an effect. 

Now the calculations of Elliott and Flowers (1955), which have enjoyed 
considerable success in this mass region, predict that the M1 transition 
from the 1:59 Mev to the 200 key levels should be near single-particle 
strength, so to produce the degree of mixing required above, the E2 
strength would have to be very much greater than that of a single particle. 
In addition their predicted value of « for the 200 kev state is 
+3-35 (+5°,) nm. (ef. the single d;), particle of the Schmidt model gives 
je=-+ 4-793 n.m.) which is in excellent agreement with our measured value 
on the assumption of no quadrupole interaction. On balance, it seems 
reasonable to infer that the simplest interpretation of our results is correct, 
that the calculations of Elliott and Flowers on the states of 1°F are still 
in accord with measurements, and that ~=-+-3-50-+0-24 n.m. is a true 
measure of the magnetic moment of the 200 kev state of 1°F. 
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Release of Stored Energy and Changes in Line Shape 
During Annealing of Deformed Nickel 


By D. MicHELL 
Division of Tribophysics, C.S.I.R.O., University of Melbourne, Australiat 


[Received December 28, 1955] 


Previous work on the broadening of the lines in the x-ray diffraction 
patterns of deformed metals has, in general, been restricted to a comparison 
of the patterns for the heavily deformed and the fully annealed states. 
However, Wilson and Thomassen (1934) investigated the x-ray line breadth 
and other properties for a number of deformed metals after annealing at 
various temperatures. For nickel, the recovery of line breadth occurred 
in two stages. Further, Clarebrough, Hargreaves and West (1955, 1956) 
have shown that the energy stored in deformed nickel is released in at least 
two stages during annealing and that there are corresponding changes in 
electrical resistivity, density and hardness. It seemed desirable to corre- 
late these changes with changes of line shape during annealing. As the 
first stage of this programme, the energy stored in deformed nickel powder 
has been measured and line shapes determined on the same powder at 
various stages of annealing. The nickel used was identical with that used 
by Clarebrough, Hargreaves and West. 

The deformed nickel powder was produced by surface grinding followed 
by sieving through a 350 mesh per inch sieve. Specimens for the measure- 
ment of stored energy and of the diffraction pattern were prepared by 
pressing this powder into suitable holders without a binder. The 
measurement of the stored energy and the manner of its release were 
determined by the calorimetric method of Clarebrough, Hargreaves, 
Michell and West (1952) using a heating rate of 6°c/min. The x-ray 
diffraction pattern between 30° and 130° in 26 was determined at room 
temperature with a Geiger counter spectrometer (monochromatic Cu Ka 
radiation) for the one specimen, as deformed and after heating at 6°c/min 
to a number of temperatures up to 870°c. 

The results of the stored energy measurements are shown in fig. 1. 
Curve (a) shows the rate of release of energy and curve (b) the energy 
remaining in the specimen, both as functions of the temperature. 
Interpreted in the manner of Clarebrough, Hargreaves and West (1955) 
for solid nickel, region A corresponds to recovery and region C to recrystal- 
lization. An intermediate region present in solid specimens is either absent 
in the filings or completely masked by the broad region C. The low 
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temperature at which the peak, region C, begins suggests that the maximum 
deformation present in the powder is greater than that in solid specimens 
of nickel deformed to fracture, while the large range of temperature over 
which this region occurs indicates that a wide range of deformations is 
present. 

Fig. 1 
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Release of energy from deformed nickel powder during annealing as a function 


of temperature. 
(a) Rate of Release of energy. (b) Energy still stored in specimen. 


The significant features of the diffraction patterns may be summarized 
as follows : 


(1) the background intensity remained constant during annealing ; 

(2) the integrated intensities of the 222 and 400 lines remained constant 
during annealing, while those of the 111 and 200 lines decreased by about 
7% over the temperature range 400-870°C ; 

(3) the tails of all diffraction lines remained practically unaltered up to 
400°c, although the peak intensities of the lines increased ; 

(4) from the measurements of line positions no evidence was found for 
stacking faults or changes in lattice parameter. 


The line shapes have been analysed by various methods (Hall 1949, 
Warren and Averbach 1950, 1952, Williamson and Smallman 1954). In 
each case the broadening has been attributed partly to strain and partly to 
particle size effects. Values of strain and ‘ apparent particle size ’ have 
been calculated as functions of temperature by each of the three methods. 
The results obtained by applying the Warren and Averbach method of 
Fourier analysis to the 111 and 222 reflections are shown in fig. 2. The 
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results for other reflections and other methods are qualitatively similar. 
Although the numerical values depend markedly on the method of analysis 
used, the decrease in strain without change in ‘ apparent particle size : 
on annealing to 400°c and the further decrease in strain with an increase in 
‘ apparent particle size ’ on annealing to 870°C are obtained by all methods 
of analysis. 


Fig. 2 


ZN 
id 
es 
_ 06 Sie ons 
— °o 
o ay, 
Da 
iS] 
6 04 1600 |N 
w 
Cc) 
Ss) 
— O 
o rer 
[= 
— 0:2 800 0 
my a 
mn 
= 
cS) 
re 
rS) 


ce) 200 400 600 
temperature (°C) 


Results of the Fourier analysis of the shape of the 111 and 222 lines from 
deformed nickel powder. 
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The energies calculated from the strains derived by the three methods 
of analysis are all considerably lower than the measured value of stored 
energy of 3-6 cal/g. The values calculated using the mean strain from the 
Warren and Averbach analysis for the 111 and 222 reflections are indicated 
in fig. 2. 

Gay, Hirsch and Kelly (1954) have suggested that a deformed metal 
consists of particles of comparatively low strain surrounded by heavily 
deformed regions. If the ‘apparent particle size’ obtained from our 
analyses can be identified with the size of their particles, the decrease in 
lattice strain without change in ‘ apparent particle size’ at temperatures 
below 400°c is consistent with their model for recovery. In addition, the 
experimental observations (1) and (2) above indicate that the diffracted 
intensity from the heavily deformed regions goes predominantly into the 
tails of the lines. Thus observation (3), which is independent of the 
analyses, is also consistent with the Gay, Hirsch and Kelly model for 
recovery. 

The stored energy measurements indicate that recrystallization 
commences at approximately 400°c. In agreement with this, as the 
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temperature of annealing is raised above 400°c, the intensity in the tails 

of the x-ray lines decreases rapidly, the calculated lattice strains decrease 

and the ‘apparent particle sizes’ increase. The decrease in integrated 

intensity of the 111 and 200 lines is probably due to an increase in extinc- 

tion resulting from an increase in both particle size and crystal perfection. 
A full account of this work will be published shortly. 
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The Spin and Magnetic Moment of 116*In 


By P. B. NurrEer 
A.S.R.E., Portsmouth, Hants. 


[Received December 9, 1955] - ~ 


A CONVENTIONAL (Davies, Nagle and Zacharias 1949) atomic-beam 
apparatus using magnetic resonance was employed to study the spinand 
magnetic moment of the 116In isomer of 52 min half-life. Transitions of 
the 4F=0, 4m=-+1 type between the magnetic sub-states of the pj/. 
atomic ground state were observed in a weak magnetic field and from them 
the nuclear spin was found to be 5. Measurements in a stronger field 
indicated a hyperfine structure splitting (Goodman and Wexler 1955) 
of Jv—8670+170 Mc/s and positive magnetic moment, the value of 
which was calculated (Fermi 1930) as 4-21 --0-08 nuclear magnetons. 
Indium metal was irradiated in the Harwell pile (BEPO); when it 
was placed in the oven it had a specific activity of about 2 curies/g mostly 
due to 116*In. To detect the beam, the atoms were ionized on the surface 
of the usual oxidized tungsten strip and then collected alternately on two 
similar targets. The radio-frequency power which caused the transitions 
was switched on and off in a regular one minute cycle as the target poten- 
tials were changed. Measuring the ratio of the activities accumulated 
on the targets after runs of about 20 minutes allowed us to determine 
resonances and eliminate slow drifts of beam intensity and source decay. 
The magnetic field was monitored continuously by observing resonances 
due to a stable Indium isotope through measuring the ion current. 
In order to compare the result with the shell model (Mayer 1950) the 
magnetic moments were calculated for possible configurations, assuming 


no configuration mixing. 
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The magnetic moments were calculated on the basis of j-j coupling 
of the odd nucleons. The value used for the magnetic moment of 
the g4. proton in these calculations was that found in the nuclei 
1130) and 145In. Experimental values for the 8/2, ds. and d;;. neutrons 
were also available, but for the g,/. and h,,/. neutrons magnetic moments 
were estimated assuming the intrinsic magnetic moment to be depressed, 
as indicated by Bloch (1951). The configurations 29/281/2, 89/245/2 and 
2 9/9 27/2 are seen to give moments in reasonable agreement with experiment. 


Configuration $9/251/2  9/2A3/2 S9/245/2 SoSr72 Sore 
Calculated 

magnetic 

moment +4-49 +5-58 +- 4:26 +4-30 + 1-37 


The author is indebted to D.S.I.R. for the provision of a maintenance 
award throughout the course of this work and also wishes to thank 
Mr. West of A.E.R.E. for his help in connection with the provision of 
active material for the experiment. 

Further thanks are extended to Dr. K. F. Smith, Mr. R. 8. Title and 
Mr. G. H. Rees for their assistance during the experiments and especially 
to the former for many previous discussions. 
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The Absorption of Sound in Liquid Helium below 1°x 


By J. A. Newewiyt and J. Winks 
Clarendon Laboratory, Oxford 


[Received March 28, 1956] 


MEASUREMENTS using a pulse technique are in progress to determine the 
absorption of sound in liquid helium as a function of pressure, and a brief 
summary of the results between 1-2°x and the lambda point has been 
given (Newell 1955). At temperatures below 1-2°x there is a large 
maximum in the absorption. and the figure shows our results in this 
region for sound of a frequency of 14-4 Me/s. (As the absorption depends 
on the density of the liquid p, and the frequency v and velocity c of the 
sound which all vary with pressure, our results are given as the coefficient 
of absorption divided by the normalizing factor 27v2/pec?.) Although 
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these results have the same general form as those obtained by Chase 
and Herlin (1954, 1955) under the vapour pressure at 12-1 Me/s, there is 
the significant difference that their maximum shows two peaks whereas 
we observe only one. These authors remark that the presence of two 
peaks gives good support to a theory of the absorption due to Khalatnikov 
(1950, 1952) ; each being associated with one of the two relaxation times 
involved in the theory. 


The coefficient of absorption of sound of 14-4 Mc/s in liquid helium at pressures 
ranging from the vapour pressure to 24-7 atmospheres. 
(As explained in the text, the coefficient of absorption 
is normalized by the factor 277v*/pc?.) 


One might attempt to resolve the discrepancy between the experi- 
mental results by recalling that on Khalatnikov’s theory the relaxation 
times vary with temperature, and therefore the position of the peaks will 
depend on the frequency of the sound. The two relaxation times vary 
with temperature at different rates, and at a certain temperature become 
equal in magnitude; thus if the absorption were measured at the 
frequency corresponding to this particular relaxation time one might 
expect to see only one peak. The difference between the frequency 
used in our experiments (14:4 Mc/s) and in those of Chase and Herlin 
(12-1 Mc/s) is too small to account for the observed difference in absorp- 
tion, but we decided that a final elucidation of this point would have to 
await the conclusion of measurements now in progress (with Dr. K. 
Dransfeld) at a frequency of 6 Mc/s. However, a recent note (Whitney 
1955) suggests that the double peak of Chase and Herlin is spurious and 
it therefore seems worthwhile to publish our data on this point. It is 
to be noted that the present results are in better accord with the theory, 
for while Chase and Herlin claim that the two peaks give evidence of two 
relaxation times, Khalatnikov’s expression for the absorption does in 
fact lead to only one peak (Arkhipov 1954, Newell 1955 b). 
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Our results may be more profitably discussed when the present measure- 
ments at different frequencies are completed and we have somewhat more 
information on the temperature dependence of the relaxation times. 
However, several points appear to be well established. (1) The peak is 
rather wide and at temperatures below that of the maximum the absorp- 
tion is much higher than the theoretical values. This may well be due to 
relaxation effects associated either with the first viscosity or with the 
scattering of phonons through small angles, neither of which processes 
have been fully considered by Khalatnikov. (2) As we will discuss at 
a later date, the variation of the height of the maximum with pressure is 
in quite good agreement with the theory, bearing in mind the approximate 
nature of some of the constants. (3) The way in which the temperature 
at which the peak occurs varies with pressure does not seem to correlate 
with the theory. 

Finally we would like to confirm the observations of Whitney regarding 
the care needed when using quartz crystals to generate pulses of sound in 
liquid helium. Helium has such a low density that the damping it 
produces: is very small, and even crystals which behave well in other 
liquids tend to ring and oscillate in spurious modes. Thus the pulse 
shape can become very irregular and the radiation pattern be so distorted 
that too large a decay is observed between successive reflections of the 
same pulse. As an example of such effects, we mention that in one 
experiment with a pair of aligned crystals acting respectively as trans- 
mitter and receiver, a very marked improvement of the pulse shape was 
obtained by inserting a diaphragm between them. Moreover, although 
the crystals were 16mm in diameter, the presence of even a 4mm 
diaphragm did not reduce the size of the signal. 
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ERRATUM 


The Creep of Cadmium Crystals at Liquid Helium Temperatures, by 
J. W. Guen, 1956, Phil. Mag., 1, 407. 


In the table the units under the headings ¢, and ¢, should, in both 
cases, read (min) not (see). 
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